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Abstract 


Thermal effects on the pin-bearing behavior of an IM7/PETI5 composite 
laminate are studied comprehensively. A hypothesis presents factors influencing 
a change in pin-bearing strength with a change in temperature for a given joint 
design. The factors include the change in the state of residual cure stress, the 
material properties, and the fastener fit with a change in temperature. 

Experiments are conducted to determine necessary lamina and laminate 
material property data for the IM7/PETI5 being utilized in this study. Lamina 
material properties are determined between the assumed stress free 
temperature of 460°F down to -200°F. Laminate strength properties are 
determined for several lay-ups at the operating temperatures of 350°F, 70°F, and 
-200°F. 

A three-dimensional finite element analysis model of a composite laminate 
subject to compressive loading is developed. Both the resin rich layer located 
between lamina and the thermal residual stresses present in the laminate due to 
curing are determined to influence the state of stress significantly. For the 
laminate modeled, the effect of modeling temperature dependent material 
properties and nonlinear stress-strain behavior was found to be negligible. 
Simply using the material properties measured at the operating temperature of 
interest was sufficient for predicting stresses accurately in a linear analysis for 
the current problem. 

Pin-bearing tests of several lay-ups were conducted to develop an 
understanding on the effect of temperature changes on the pin-bearing behavior 
of the material. As expected, for all lay-ups that failed in the bearing mode, pin- 
bearing strength decreased with an increase in temperature. Micrographs of 
failed specimens revealed severe damage in the outermost two plies with shear 
cracks emanating from that region. 

A computational study investigating the factors influencing pin-bearing 
strength was performed. A finite element model was developed and used to 
determine the residual thermal cure stresses in the laminate containing a hole. 
Very high interlaminar stress concentrations were observed two elements away 
from the hole boundary at all three operating temperatures. The interlaminar 
stresses quickly reduced to zero within one laminate thickness from the hole 
boundary. Changes in fastener fit with a change in temperature were predicted 
to be insignificant on influencing changes in pin-bearing strength. 

The pin-bearing problem was modeled assuming a rigid frictionless pin 
and restraining only radial displacements at the hole boundary. A uniform 
negative pressure load was then applied to the straight end of the model. A 
solution, where thermal residual stresses were combined with the state of stress 
due to pin-bearing loads was evaluated. The presence of thermal residual 
stresses intensified the interlaminar stresses predicted at the hole boundary in 
the pin-bearing problem. This dissertation shows that changes in material 
properties drives pin-bearing strength degradation with increasing temperature. 
The thermal residual cure stresses affect the extent that the pin-bearing strength 



changes with temperature, where they can only lessen the strength degradation 
with increasing temperature. 
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Chapter 1 
Introduction 

The desire for lighter weight aircraft and spacecraft has led to the 
utilization of polymer matrix composite materials for many structural components 
in advanced vehicle designs. Laminated composite panels have the ability to 
carry high mechanical loads while remaining lightweight in comparison to metallic 
panels. In advanced vehicle structural applications, the composite panels may 
be exposed to temperature extremes. This dissertation considers a critical area 
for application of composites at temperature extremes: the problem of 

mechanically fastened joints. 

1.1 Background 

To provide the rationale for the dissertation research, the necessity of 
mechanically fastened composite joints is briefly discussed. Next, the loading of 
composite joints at elevated and cryogenic temperatures is described followed by 
the topic of pin-bearing loads on a composite laminate and associated failure 
modes. 

1.1.1 Mechanically Fastened Composite Joints 

Size limitations in large-scale structures fabrication in addition to economic 
factors require that structural components be fabricated as components with 
subsequent assembly to produce a final vehicle. Due to the inability of adhesive 
bonds to transfer high mechanical loads between structural components, 
mechanically fastened joints are often used in vehicle design. A schematic 
drawing of a mechanically fastened composite joint concept as a component in 
an aircraft wing is illustrated in Figure 1.1. The joint concept shows two 
honeycomb sandwich panels with composite facesheets joined to a titanium 
splice plate through two rows of fasteners on each panel. By providing a load 
path through bolts and splice plates, the bolted splice joints transfer mechanical 
loads between structural components. 

The need to join composite panels through bolted splice joints has raised 
many issues not encountered with metallic joints. Anisotropic, nonhomogeneous 
and viscoelastic properties of composite laminates in addition to their brittle 
nature require special attention in developing structurally reliable and efficient 
joints. The joining of composite laminates through bolted splice joints 
consequently became a major topic of research from the 1970s to the present. 
In 1986, the first U.S. government design guide was published that provides 
analysis tools (computer codes) and guidance for the design of composite bolted 
joints. 11 


1.1.2 Thermo-mechanical Loading 

When a structure is exposed to a change in operating temperature, 
thermally induced stresses arise from constrained thermal expansion or 
contraction. In addition, both thermal and mechanical properties may change 
over the range of temperature. The combined effect of altered thermal 
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environment and additionally applied mechanical loading results in a complex 
thermo-mechanical loading of the structure. 

Elevated Temperatures. During typical aircraft flight maneuvers, bending loads 
are induced along the span of the wing. Wing bending results in significant in- 
plane axial loads for the upper and lower wing surface panels. The loads may 
either be tension or compression; for example, when the upper surface is in 
compression, the lower surface is in tension. The axial load is depicted in Figure 
1.2 for a surface wing panel joint experiencing a tensile load. A side view 
through the thickness of the wing joint is displayed in the figure with the applied 
mechanical load, N x . For high speed aircraft structures, like that shown in Figure 
1.1, the outer surface structural joints are also exposed to aerodynamic heating. 
For a Mach 2.5 civil transport, temperatures have been shown to reach 350 F on 
the outer surface of the aircraft. 1 2 Thermal effects in addition to the mechanically 
applied loads create a state of thermo-mechanical loading for the surface splice 
joint. The elevated temperatures may degrade the load carrying capabilities of 
the joint and consequently impose additional challenges in the joint design. 
Unfortunately, the design tools of Ramkumar et al. do not include temperature 
effects. 1 1 

Cryogenic Temperatures. On the other end of the spectrum, there is the 
possibility that laminated composite bolted splice joints will be necessary in 
cryogenic environments. Research studies have shown the necessity for 
composites as primary structure for a cost effective reusable launch vehicle(RLV) 
to replace the Space Shuttle. 1 3l 14 A schematic of a composite bolted splice joint 
in a RLV concept is illustrated in Figure 1 .3. The figure shows an RLV concept 
with internal liquid oxygen (L02) and liquid hydrogen (LH2) cryogenic tanks and 
an intertank structure located between the cryogenic tanks. In such a design, 
a composite bolted splice joint is used as an assembly joint in mating the 
intertank structure with the cryogenic tank structures. As shown in the figure, the 
joint is subject to axial loading, N x . Due to the joints proximity to the cryogenic 
tanks, the joint will be exposed to cryogenic temperatures. Consequently, the 
effect of cryogenic temperatures on the joint behavior requires consideration in 
the joint design. 

1.1.3 Composite Laminate Subject to Pin-Bearing Loading 

When the bolted splice joints shown in Figures 1 .2 and 1 .3 are subject to 
in-plane mechanical loads, the composite laminate facesheets experience 
complex loading, where bearing loads are a significant component. The bearing 
load is compressive due to the fastener bearing on a finite area of the composite 
laminate at the bolt hole boundary. 

A simplistic case of pin-bearing where there is no out-of-plane compression of 
the composite laminate due to tightening of the bolt is shown in Figure 1 .4. The 
composite laminate bearing specimen has width, w, edge distance, e, bolt hole 
diameter, D, and thickness, t. The figure also shows a tensile load, P, being 
transferred to the composite laminate through the pin with diameter, d. Due to 
the load transfer through the pin-laminate contact, the laminate is subject to 
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compressive bearing stresses in the vicinity of the bearing surface. The contact 
area varies as the load P increases. 

A typical load-deformation curve obtained from a bearing test is shown in 
Figure 1.5. 15 For bearing testing, the deformation, 5, is measured across the 
bolt hole and is referred to as the bolt hole elongation. On the curve, the 
proportional limit represents the onset of nonlinear material behavior. The yield 
strength corresponds to 4% hole elongation in the composite laminate, and the 
ultimate load represents the maximum sustained load of the joint prior to failure. 
The bearing strength F br of the laminate is calculated by: 


F br = P/tD (1.1) 

The subscripts bry and bru are also used to differentiate between yield and 
ultimate bearing strength, respectively. 

The four basic laminate failure modes associated with composite bolted 
joints are illustrated in Figure I.6. 1 5 The tension failure mode is characterized by 
a fracture of the laminate across its width from the hole to the edges, the 
shearout failure mode is characterized by a “pull-out” fracture between the hole 
and laminate end, and the cleavage-tension failure mode is characterized by a 
simultaneous fracture across the width to one edge and between the hole and 
the laminate end. The bearing failure mode is observed as a localized crushing 
of the laminate in front of the bearing surface. The type of failure mode observed 
during bearing testing is affected by the w/D and e/D ratios of the bearing 
specimens, which can easily be explained through strength of materials 
considerations. For a given laminate, the tension failure mode is associated with 
insufficient w/D ratios, where the specimen is not wide enough to prevent a 
tension failure. The shearout failure mode is dependent on the e/D ratio, where e 
can be made larger to prevent shearout failure. The cleavage-tension failure 
mode is a combined tension and shearout failure when the tension strength 
corresponds to the shear strength of a specimen. For sufficiently large e/D and 
w/D ratios, the bearing failure mode will result. 

1.2 Previous Research 

Previous research pertinent to this dissertation is reviewed in this section. 
Research on thermal effects on bearing behavior of composite joints is presented 
by first reviewing experimental investigations, followed by analytical studies. 
Residual stress concentrations at bolt holes is then reviewed. Finally, issues 
concerning previous studies are revealed. 

1 .2 .1 Experimental Investigations on Bearing Behavior 

Previous experimental investigations on pin-bearing strength, bolt bearing 
strength, and the effect of clamp-up forces for composite laminates with varying 
operating temperatures are reviewed. 

Pin-bearing Strength. Three experimental investigations reported on the effect of 
temperature on the pin-bearing strength and failure mode for particular 
composite laminates. Pin bearing tests were conducted by Kim et al. (1976) at 
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room temperature and at 260°F on three different graphite/epoxy lay-ups. 1 6 
About a 30 percent reduction in the pin bearing strength was observed for a 
temperature increase to 260°F for all three lay-ups. No change in the bearing 
failure mode was observed at the elevated temperature. Another study was 
conducted by Wilkins in 1977 which also investigated the effect of temperature 
on joint strength for the same graphite/epoxy material for five different lay-ups. 
For all specimens, a moderate temperature increase from room temperature to 
200°F showed no significant effect on the specimen strength. However, unlike 
Kim et al.’s study where the bearing specimens failed in a bearing mode, the 
bearing specimens here failed in a net tension mode. In contrast, a pin bearing 
strength reduction of approximately 40% with a temperature increase to 350 F 
from room temperature was reported by Wilson in 1979 for graphite/epoxy 
specimens with a range of geometric parameters. 1 8 Interestingly, for the 0 
degree dominated lay-up tested, where the strength difference was pronounced 
for varying specimen width at room temperature, the strength became 
independent of specimen width at 350°F. 

Bolt Bearing Strength. Several references encompassing both elevated 
temperature and cryogenic strength testing of graphite/polyimide composite 
bolted splice joint specimens were found in the literature. In 1979, Perry and 
Hyer established the testing procedures in their study. 1 Testing of 16-ply double 
lap quasi-isotropic specimens were conducted at temperatures of -250 F, 75 F, 
and 600°F. Wichorek observed a 30% reduction in bearing strength at the 600 F 
testing temperature in comparison to room temperature, and an 18% increase in 
bearing strength at the low temperature of -250°F in comparison to room 
temperature. 1 10,1 11 Other experimental investigations also revealed strength 
reductions at elevated temperatures for bolted composite material joints. 
i.i4, i.i5, i.i6, i.i7 0ne study concluded that fracture load and mode can be 

predicted for joint geometry and temperature for a specific material system, 
where the authors proposed an empirical equation to estimate fracture load.^ 
However, for carbon fiber composite joints tested at temperatures up to 285 F, 
negligible strength reductions were observed. 1 19 In general, past research 
suggests that the magnitude of strength reduction is at least somewhat 
dependent on the composite unnotched material properties variations with 

temperature. , . 

The effect of lay-up on bearing strength at room temperature (RTD) and at 

a 250°F elevated temperature wet condition (ETW) was investigated by Garbo 
and Ogonowski for graphite-epoxy bolted joints. 1 20 Three lay-ups of varying 
percentage of 45° plies were tested in both tension and compression for a double 
shear two-fastener in-tandem specimen with a 50 in-lb torque. A 30-40 ksi loss in 
ultimate bearing strength was observed for all three lay-ups tested at the ETW 
condition. This result is in contrast to other reported pin-bearing strength results 
where the magnitude of the strength reduction was dependent on the lay-up of 
the material being tested. 1 6 ' 17, 1 Additionally, the effect of edge distance, e, 
and specimen width, w, (see Figure 1.4) on the bearing strength was also 
investigated by Garbo and Ogonowski. For a pin diameter d, they observed an 
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insensitivity to both e/d and w/d ratios, with a reduced bearing strength capacity 
for the ETW condition. 

An extensive test program was conducted with the objective of providing 
information on the bearing behavior of graphite cloth epoxy in addition to quasi- 
isotropic tape laminate bolted joints at temperatures up to 450K. 1 ' 21, 1 22 Strength 
degradations at elevated temperatures were observed for single and 
multifastener joint tests. The empirical approaches of Hart-Smith were utilized to 
explain elevated temperature test results. The techniques of Hart-Smith were 
developed from room temperature testing where stress concentration factors are 
computed based on composite bolted joint strength data from joints failing in 
tension. 1 23 1 24 Baile et al.’s test results from tension failure specimens showed 
that raising the temperature produced significant increases in the stress 
concentration. This was explained to be due to higher temperatures reducing the 
resin strength and stiffness, thus lessening it’s ability to support the fibers and 
transfer the load adjacent to the hole. Additionally, when tension failures 
occurred, there was little if any decrease in the joint stiffness as a function of 
temperature, but when bearing failures occurred, there was a clear decrease in 
joint stiffness over the linear load-deflection range. 

The mechanical behavior of quasi-isotropic bolted joints was examined 
and discussed by Meakawa et al. in 1985, where test temperature, specimen 
width, and edge distance were all varied with zero bolt-torque clamping. 125 
Examination of the results shows that at 120°C, all failure modes changed to 
bearing with a lower bearing strength than at room temperature. The change to a 
bearing failure mode was attributed to the actual decrease in the compressive 
strength for the laminate with increasing temperature. The width and edge 
distances required to raise the joint strength at room temperature did not provide 
any significant strength improvement at elevated temperatures where bearing 
failures occurred. In retrospect, these results are similar to those observed by 
Garbo and Ogonowski for graphite/epoxy bolted joints, where increasing 
specimen width and edge distance did not provide the increased strength 
capacity for the ETW case as it did for the RTD case. 1 20 Examination of the 
material properties of the lamina reported by Garbo and Ogonowski shows a 
53% reduction in the fiber direction lamina compressive strength and a 33% 
reduction in the matrix direction lamina compression strength at the ETW 
condition and no change in the fiber direction tension ultimate strength with no 
reporting on the matrix direction tension strength at the ETW condition. 

Bolt Clamp-up Effects. The beneficial effects of bolt clamp-up forces in 
increasing joint strength has been well documented for room temperature 
composite joints. 1 26, 1 7 ' 1 28 The increase in strength can be attributed to a 
combined effect of increasing the portion of load being transferred by friction at 
the joint, thus decreasing the bearing load at the bolt hole and also the lateral 
constraint increasing bearing strength by reducing out-of-plane displacement 
under contact and hence reducing the out-of-plane shear stress. 1 29 The effect of 
temperature on joint strength for varying clamp-up force was investigated by 
Horn and Schmitt in 1994. 1 0 They conducted tests of single-shear thermoplastic 
composite specimens with the objective of establishing a correlation between 
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bearing strength and clamp-up force. Their test results showed that increasing 
temperature decreased the bearing strength proportionally over the range of 
clamp-up force. Both the room temperature and elevated temperature results 
showed the same trend of increasing bearing strength with increasing clamp-up 
force. Wright et al.(1997) investigated the time- dependent bolt bearing behavior 
of thermoplastic joints at elevated temperature. 1 31 They observed time- 
dependent behavior of IM7/K3B joints only in a narrow band of loadings for an 
unclamped condition. For the clamp-up tests, no bearing creep behavior was 
observed up to 1 50% of the baseline load for the joint. Clamp-up force appeared 
to have the most significant effect on coupon bearing strength at temperature. 
They also conducted bearing bypass interaction static tests at temperature. The 
test results showed more significant reductions in strength under compressive 
loading at the elevated temperature than under tensile loading conditions. 

Conclusions. The experimental investigations presented show that an increase 
in temperature leads to bearing strength reductions in many polymer matrix 
composite laminates. The extent of the bearing strength reduction has been 
shown to be affected by not only the temperature change, but also by the 
material system and lay-up being evaluated. The extent of bearing strength 
reduction with increasing temperature was also shown to be independent of 
clamp-up force for two material systems. Additionally researchers observed that 
the bearing failure mode is more likely to occur than the tension failure mode at 
elevated temperature due to a larger reduction in the laminate compressive 
strength at the elevated temperature in comparison to the laminate tensile 
strength. Consequently, based on the elevated temperature laminate tensile 
strength and bearing strength data for the bearing failure mode, one could easily 
design a joint with an optimum width, or distance between fasteners, for elevated 
temperature applications. 

Additionally, a study was presented that showed increase bearing strength 
at the low test temperature of -250°F for graphite/polyimide specimens. 
Unfortunately, the material properties of the quasi-isotropic laminate tested were 
only reported at room temperature. 

1 .2.2 Analytical Studies on Bearing Behavior 

Finite element analysis methods have been applied to investigate the 
bearing behavior of composite joints at elevated temperatures. A two- 

dimensional, in-plane, finite element analysis was employed by Maekawa et al. in 
1985 to determine principal stresses in the vicinity of the bolt hole under bearing 
loading conditions. 1 25 By examining risk factors computed from the ratio of 
principal stresses to the compressive and tensile laminate strengths, the ability to 
predict the failure mode proved to be feasible. As the temperature increased, 
higher compressive risk factors were observed at the bolt hole, thus increasing 
the risk for a compressive failure mode. The risk assessment compared well with 
test results where bearing failures were more prominent at elevated 
temperatures. 

An incremental finite element code was used by Chen and Lee to study 
the effect of temperature on the contact stresses at pin loaded holes in 
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composite laminates. 1 32, 1 33 A two-dimensional finite element plate model of the 
composite laminate with a constrained pin, as displayed in Figure 1.7, was 
developed for their investigations. They adopted a quadratic failure criterion to 
predict failure and study the failure process at different temperatures and for 
different lay-ups and materials. Their analysis results showed bearing strength 
degradations with increasing temperature, and that the extent of strength 
degradation with increasing temperature was dependent on the material system 
and lay-up being analyzed. Once again, analysis results showed a change in the 
failure mode with temperature, from that of a tension failure mode at room 
temperature to a bearing failure mode at elevated temperatures. In Chen and 
Lee, their analysis model was first subject to a temperature increase while a 
simply-supported boundary condition was imposed on the boundary where the 
mechanical load, U e , was subsequently applied. 133 Initial thermally induced 
contact loads were observed at the hole boundary with the temperature increase 
prior to the introduction of the mechanical load. Increasing the pin clearance 
resulted in a decrease in the initial thermal load as would be expected for their 
joint model. The bearing capacity at temperature due to mechanical loading, 
which they defined as the residual bearing strength, consequently was strongly 
degraded by the initial thermal load. 

Conclusions. In the literature presented, two-dimensional finite element analyses 
have been used in evaluating bearing behavior. The analysis results provide 
insight on the effect of temperature on the bearing failure mode, which was in 
agreement with experimental observations on the increased likelihood of a 
compressive failure mode at elevated temperatures. Analysis results also 
confirmed the dependency of strength degradation with increasing temperature 
on the material and lay-up being considered. 

1.2.3 Residual Stress Concentrations at Bolt Holes 

During the fabrication of a composite laminate, when layers are oriented in 
varying fiber directions, thermal residual stresses arise due to curing the 
orthotropic layers at high elevated temperatures and then cooling in a bonded 
state. The existence of the residual thermal stresses in composite laminates 
leads to thermal residual stress concentrations when bolt holes (or circular 
cutouts) are introduced into the laminate. 

Three papers were found in the literature that included thermal residual 
stresses at bolt holes in composite laminates in their investigations. In 1976, 
Rybicki and Schmueser performed a three-dimensional finite element analysis on 
a composite laminate containing a circular hole. 134 A temperature loading 
condition of -200°F was applied for several lay-ups of graphite/epoxy laminates. 
Constant material properties were assumed for their linear analysis. Both 
significant in-plane and tensile transverse stress results were presented at the 
bolt hole only. No comparison was made to the stresses that would be present in 
the absence of a bolt hole. A nonlinear finite element approach was used by 
Chang in 1990 to study the stress concentrations at cutouts in asymmetric 
laminates. 1 35 His two-dimensional plate analysis model included temperature- 
dependent material properties and accounted for nonlinear geometric effects. 
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Large tensile tangential stress concentrations were observed for the various lay- 
ups studied. Stress concentrations were also shown to increase with decreasing 
hole diameter. In 1986, Smith and Pascal included the contribution of thermal 
residual stresses in estimating the interlaminar normal stress at a bolt hole 
subject to a pin-bearing load. 136 An approximate expression for interlaminar 
normal stress due to bearing loads was derived based on equilibrium 
considerations. Their simple estimate of the total interlaminar normal stress at a 
bolt hole showed a significant contribution due to thermal residual stresses. 

These papers identify, through analysis, the existence of thermal residual 
stress concentrations in the vicinity of bolt holes in composite laminates. 
Consequently, these residual stress concentrations may affect the bearing 
strength of a composite laminate. Unfortunately, to this author s knowledge, no 
attempt has been made to measure the residual stress concentrations at the bolt 
hole in a composite laminate. This lack of data is due to the inability to measure 
strain accurately in the vicinity of very localized large gradients. 

12 4 Issugs 

Although - the literature has presented reduced joint strengths at elevated 
temperature conditions for many polymer matrix composite material systems, the 
actual mechanisms attributing to the bearing strength degradation have not been 
well defined. Additionally, the improved bearing strength that was observed at a 
low temperature requires further investigation as to the nature of the strength 
change. A basic understanding of thermal effects in a composite joint is key to 
producing a structurally efficient and reliable composite joint for high and low 
temperature operating environments. Knowledge of the factors attributing to the 
variation in strength can be used to improve joint designs for extreme 
temperature applications. 

The literature has identified the possibility of significant thermal residual 
stresses at bolt holes in composite laminates. Although the residual thermal 
stresses were shown to have a significant contribution to the interlaminar normal 
stress at a bolt hole subject to bearing loading in one very simplified 
investigation 1 36 , to this author’s knowledge, the contribution the thermal residual 
stresses have in affecting the bearing strength has not been investigated. 
Furthermore, the existence of thermal residual stress concentrations has not 
been considered in any investigation found in the literature concerning the effect 
of operating temperature changes on bearing strength. A detailed three- 
dimensional analysis on the state of residual stress due to curing and the effect 
of temperature change coupled with bearing loads is necessary to provide more 
insight on the effect of temperature on bearing strength of composites. Also, with 
the recent development of materials cured at higher temperatures, thermal 
residual stresses become more significant and consequently require 
consideration in evaluating both laminate and bearing strength changes with 
changing temperature. 
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1.3 Hypothesis 

When a composite laminate containing a bolted joint is subject to a 
uniform change in temperature, AT, the author hypothesizes that the pin-bearing 
strength change, AF br (AT), is due to the combined effect of three basic changes 
in the joint? -Tire three basic changes are: 

1 . a change, AS, in the state of residual thermal stress at the hole boundary in 
the composite laminate due to a change in temperature, 

2. a change, AM, in the material properties of the constituent materials, and 

3. a change, AB, in the fastener fit due to changes in the bolt hole geometry and 
radial expansion or contraction of the fastener. 

Thus, an expression for the change in pin-bearing strength as a function of a 
change in temperature is proposed as 

AF br (AT) = f (AS ,AM, AB) (1.2) 

where AS ,AM, and AB are all dependent on AT for a given joint design. The 
significance each change, AS , AM, and AB, has on influencing the bearing 
strength may depend on the material system and joint configuration being 
considered. Consequently, the expression may need to be revisited when any 
change in a joint design occurs to determine what the significant factors are for 
that particular joint. 


1 .4 Goal and Objectives 

The goal of the research is to develop a basic understanding of thermal 
effects on the pin-bearing behavior of a composite laminate. By examining in 
detail the extent the hypothesized changes of residual stress, material properties 
and fastener fit have on influencing the bearing behavior, a fundamental 
knowledge on the nature of the bearing strength change will be achieved. This 
knowledge can consequently be used to design more efficient and reliable flight 
vehicle joints. 

The objectives of this dissertation are to: 

1) characterize a material system within the applicable temperature range with 
the detail required to model structural response accurately under thermo- 
mechanical loading, 

2) conduct pin-bearing experiments to provide data for evaluating and developing 
a general pin-bearing analysis model for analysis under thermo-mechanical 
loads, and 

3) develop computational models with adequate detail to study the significance 
the hypothesized changes have in influencing the state of stress and 
consequently bearing behavior for the pin-bearing configuration. 


1.5 Scope 

The composite material system characterization is described in Chapter 2. 
The analysis of the composite laminate to establish laminate attributes necessary 
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in the pin-bearing finite element analyses is presented in Chapter 3. The pin 
bearing test program for validation of the computational models is described in 
Chapter 4. The computational models for investigating thermal effects on stress 
concentrations at bolt holes, fastener fit, and the state of stress under pin-bearing 
loading are described in Chapter 5. Finally, the completed research is discussed 
and conclusions are made in Chapter 6. 
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Figure 1.1. A mechanically fastened composite joint concept as a component in 
an aircraft wing. 



— 

— 



z* 

zx: 

n 


— 

— 

— 

1 

1 

1 

1 

1 

1 

i 


1 

i 

1 

i 

1 

i 

1 

1 


T5 

nnr 

r o iut 1 



Figure 1 .2. A composite wing panel splice joint concept subject to tensile axial 
loading. 
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Figure 1.3. A mechanically fastened composite joint concept as a component in 
a reusable launch vehicles cryogenic intertank structure. 
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Figure 1 .4. Pin-bearing tension test configuration. 


17 


Load, P 



0.04D 


Deformation, 8 



Figure 1 .5. Typical load-displacement curve for a bearing test. 
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Figure 1.6. Four possible composite laminate failure modes associated with 
bolted joints. 
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ite element model of plate and pin, Ref. 1.33. 
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Chapter 2 

Material System Characterization 

In this chapter, a composite laminate is first described, followed by a 
description of the IM7/PETI5 material system being utilized for this study. The 
material properties for characterization of the material system are defined, which 
includes the temperature-dependent lamina properties, the nonlinear stress- 
strain behavior of the lamina, and finally the laminate strength properties. 


2.1 Composite Laminate Description 

A polymer matrix composite (PMC) laminate is composed of layers 
referred to as lamina. Each lamina consists of high strength fibers embedded in 
a lightweight polymer matrix resin material as illustrated in Figure 2.1. Also 
shown in the figure are the principal material directions. The 1 -direction is 
parallel to the unidirectional fibers, and is also considered the fiber direction, 
since the properties in the 1 -direction are dominated by the high strength fibers. 
The 2-direction is considered the matrix direction and is perpendicular to the fiber 
direction and in-plane to the lamina. The 3-direction is also perpendicular to the 
fiber and is out-of-plane to the lamina. The lamina material is subsequently used 
to construct laminates, where the layers of lamina are stacked in a lay-up with 
various orientations of the principal material directions. The direction of the 
layers and the stacking sequence are chosen to tailor the mechanical properties 
of the composite laminate. An example of an 8-ply laminate with a (+45, 
-45 0 90) s lay-up where the subscript s denotes symmetric about the mid-plane 
z =o, is illustrated in Figure 2.2. Only the top half, or the top four plies, of the 8- 
ply laminate is shown in the figure. The laminate x, y, and z coordinate system is 
shown in the figure along with the orientations of the top four layers. Once 
stacked, the lay-up is cured in an autoclave at elevated temperatures and 
pressures to form a permanent bond between layers. 

2.2 IM7/PETI5 Material System # 

The composite material being utilized for this study is a high strength PMC 
material designated IM7/PETI5 21 . The material consists of IM7 graphite fibers 
with a PETI5 polyimide matrix material. A 100 lb batch of the lamina material in 
prepreg tape form was purchased from Cytec Fiberite for the study. 

The IM7/PETI5 material system was chosen since the system offers 
potential for use in high-speed aircraft structural applications up to temperatures 
of 350°F. Additionally, IM7/PETI5 laminates are cured at 700°F, a temperature 
much greater then the approximately 300°F temperature used to cure most 
graphite/epoxy laminates. Due to the large temperature change in processing 
the laminates, there is potential for significant thermal residual stresses due to 
curing to be present in the laminates at their operating temperatures. The 
magnitude of the thermal cure stresses in IM7/PETI5 laminates has not been 
previously investigated and is consequently being investigated here to provide 
the necessary input to evaluating their effect on bearing behavior. 
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2.3 IM7/PETI5 Material Properties 

Due to the lack of all necessary property data and the need for accurate 
material properties for the computational study, relevant material property data 
has been determined from the current batch of IM7/PETI5 prepreg. Both lamina 
properties and properties for the laminates of interest for the bearing study have 
been measured. The lamina properties are used in the analysis models for 
predicting the state of stress in the composite laminate. To aid in the 
development of the computational models, the laminate strength properties were 
measured for qualitative comparisons between lamina and laminate strength. 
Details on the property tests conducted herein are given in the Appendix. 

2,3.1 Lamina Properties 

The temperature-dependent orthotropic lamina properties obtained from 
testing in the IM7/PETI5 material are given in Table 2 . 1 . The properties are given 
at the three operating temperatures of -200°F, 70°F, and 350°F and an upper 
bound temperature of 495°F. Since the variation in E\ i c obtained from the tests 
at the different temperatures was negligible, the average £ 11c from all test 
temperatures is presented in the table and it is concluded that £ 11c is 
independent of temperature. For characterization of the material for non-linear 
analysis purposes including temperature dependent properties, the properties 
are assumed to vary with temperature piecewise linearly between the four given 
temperature values. The glass transition temperature, T g , of the material was 
measured to be T g = 460°F, after the properties were already determined up to 
495°F. 

For predicting residual cure stresses in the laminates under investigation 
in this dissertation, the T g is assumed to be the stress free temperature(SFT) 
below which thermal curing stresses develop in the laminates. The choice of the 
T g as the SFT is in accord with previous investigations on curing stresses in 
graphite/epoxv laminates, where a SFT of 350°F, the typical T g of graphite/epoxy, 
was chosen . 22 23 Therefore material stiffness and coefficient of thermal 
expansion properties are only utilized up to the T g of 460°F, for predicting 
residual cure stresses in the laminates of interest. 

The additional out-of-plane properties needed for a three-dimensional 
characterization of the composite material include E 33 , G 3U G 23 , v 31 , and a 3 . 
These properties are assumed related to the in-plane properties as follows, 

e 33 — E22 

G 31 = G 12 

F 

r - ^22 

2(1 + v 2} ) 

E „ 

V 31 — — V 12 where vi 3 = vi 2 

E n 

a 3 = a 2 

and are determined from the values given in Table 2.1. 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 
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The properties of the resin rich layer between lamina, which will be 
discussed later in Chapter 3, are assumed to resemble the m atrix dir ®c™ 
properties of the IM7/PETI5 material, which are dominated by the PETI 5 
properties. The PETI5 resin layer is hence assumed to be an isotropic material 
with the following two independent properties, 


£ r = £ 22 
v r = 0.4 


( 2 . 6 ) 

(2.7) 


where for isotropic materials, 


The simplifying assumption of constant material properties allows for a 
linear static analysis to be performed. When the assumption of constant 
properties is used to determine the stresses at a particular operating temperature 
within this dissertation, the properties at the operating temperatures given in 
Table 2.1 are all used along with the properties determined above at the 
operating temperature of interest. 


2.3.2 Non-linear Stress-Strain Behavior 

The moduli presented in Table 2.1 have been determined as a constant 
from a linear approximation of the slope of the stress-strain curve in accord with 
the procedures of the ASTM standard (see Appendix). However, the test data 
revealed that at all test temperature, G 12 , and at the test temperature of 495 F 
and 350°F, £ 220 . were not constant and alternatively were a function of strain. 
Using fourth order polynomial curve fitting of the actual stress-strain test data 
from the replicate specimens, an average curve was computed for evaluation. 
The average stress-strain curves derived from the tests conducted to determine 
G 12 and E 220 are plotted in Figures 2.3 and 2.4, respectively. The tangent moduli 
computed from these curves at specific values of strain are presented in Table 
2 2 and Table 2.3. Furthermore, the properties £33, G 23 , G31, V31, E t and G r , being 
directly related to either G 12 or £ 22 through Equations 2.1 -2.7, are also strain 
dependent. Tables 2.2 and 2.3 along with Equations 2.1-2.7 characterize the 
material’s nonlinear stress-strain behavior. The same characterization is 
assumed for both positive and negative values of shear strain. 

As can be observed from Figures 2.3 and 2.4, the most significant 
nonlinear behavior is seen at the upper bound temperature of 495°F and at only 
relatively large stresses at the other temperatures. Since during curing, the 
magnitude of stress is assumed to increases from zero at the stress-free 
temperature of 460°F, and since the nonlinear stress-strain behavior dominates 
at only large stresses and high temperatures, no analyses with nonlinear stress- 
strain properties will be pursued in the context of this dissertation. Additional 
justification on the neglect of nonlinear material behavior will be discussed later 
in Chapter 3, in light of the analysis results. All finite element analyses within the 
context of the dissertation will assume linear elastic material properties. 
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2.3.3 Laminate Strength Properties 

The measured laminate compressive strength, Fxx C , of the particular 
laminates of interest is presented in Table 2.4. The lay-ups of the laminates 
correspond to the lay-ups that will be evaluated in the bearing study. The 
compressive strength data will be used for qualitative comparisons with analysis 
result. 
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Table 2.1 Experimental Data-Lamina Properties. 


Material 

Property 

Test Temperature (F) 


-200 

70 

350 

495 

E nc(Msi) 

20.4 

20.4 

20.4 

20.4 

F\ ic(ksi) 

238. 

225. 

185. 


E i ]y(Msi) 

21.9* 

21.9 

22.9 

25.8 

Fuj(ksi) 

336. 

307. 

307. 


Enc (Msi) 

2. 

2.2 

1.5 

1.2 

Fncfai) 

47. 

37. 

26.8 


E22T (Msi) 

1.6 

1.4 

.97 

.6 

E 22T (ksi) 

6.6 1 

9. 

7.3 


G\i(Msi) 

1.04 

.92 

.69 

.62 

F\i(ksi) 

17.8 

14.6 

6.6 


V12 

.34* 

.34 

.33 

.34 

V23 

.4 

.64 

.52 

.42 

/'’iz(ksi) - ILSS 

21.3 

17.3 

12. 


Fzz(ksi) - ILNS 

5.2 

4.5 

2.7 


ai(pin/in/F) 

-.24 

-.09 

.9 

.83 + 

a2(|iin/in/F) 

10.1 

15.5 

16.7 

17.3 + 


data unavailable, extrapolated from closest value 

value given is at 460°F (a changed from value given to 0. at T=460°F) 








Table 2.4 Experimental Data - Laminate Compressive Strength Properties 


LAY-UP 

Test 

Temperature 

(°F) 

Laminate ! 
Strength 

Txxc 

(ksi) 

1 

(90,0) 8s 

-200 

144. 

70 

126. 

350 

95. 

3 

(+45,-45,0,90)4s 

-200 

118. 

70 

105. 

350 

76. 

4 

(+45,0,-45,90)4s 

-200 

115. 

70 

93. 

350 

69. 

5 

(+45, 02,-45,02, +45, 0 2 , -45, 0 2 ,+45,90 2 ,-45) s 

-200 

144. 

70 

126. 

350 

79. 



Figure 2.1 . Unidirectional fiber reinforced lamina. 
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Figure 2.2. Composite laminate construction. 







T= 350°F 



Figure 2.4. Nonlinear Stress-Strain curves for compressive matrix direction 
modulus, E 2 2 c. 
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Chapter 3 

Analysis of Composite Laminates 

An understanding of the structural behavior of the unnotched composite 
laminate is essential to investigate pin-bearing behavior. The details required in 
the model of a composite laminate to predict behavior accurately need to be 
incorporated in the pin-bearing model to understand the bearing response. 

In the following sections, the stresses in composite laminates, as pertains 
to this dissertation, are derived using classical lamination theory. The concepts 
of interlaminar stresses and residual thermal curing stresses follow. A verification 
model is then presented followed by the details of the three-dimensional analysis 
model development of an IM7/PETI5 composite laminate. 

3.1 Classical Lamination Theory 

Using the lamina orthotropic material properties, laminate behavior can 
be predicted using Classical Lamination Theory (CLT) and existing failure 
theories. 3 1 CLT assumes the lamina to be in a state of plane stress where the 
out of plane stresses, o 3 , Xi 3 and x 23 , are assumed to be zero. For an orthotropic 
lamina, the thermoelastic stress-strain relations in principal material coordinates 
are 


V 


0., 0.2 0 ' 


e, - a, AT 

ct 2 

► = 

012 022 0 


e 2 -a 2 &T • 

?\1 . 


. 0 0 0 M 


Yn 


or, equivalently, M l2 = [0X{«} 12 “Wn Ar ) 

where the 0j’s are the reduced stiffnesses for the plane stress state. 

the engineering constants, 


0 , .= 


En 


l-Wii 


022 ~ 


_ 22 


In terms of 


(3. 2) 


0.2 = 


^12 E22 


-V' 12 V 2I 


066 ' f -^12 


Figure 3.1 illustrates the transformation between the material and laminate 
coordinate systems. The transformation matrix, [T] is written as 

cos 2 # sin 2 # 2sin#cos# 

[t]=\ sin 2 # cos 2 # -2sin#cos# | (3.3) 

-sin#cos# sin#cos# cos 2 #-sin 2 # 


When relation (3.1) is transformed into the laminate coordinate system, the 
lamima stresses in the kth ply of the laminate are 
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QjJ_ Qn_ Qi6_ E x~ a , AT 

Qn_ Q^_ K-Vrl (3.4) 

Ql» Ql6 Q 6 6 y*y — a xy AT 

or, equivalently, M* = 0 ({«}* - AT? 

where \Q ] is the transformed reduced stiffness matrix and is related to [ Q ] by 

IQ] = [T] 1 [Q)[T}- J (3.5) 

and (3.6) 

CLT assumes that the laminate is thin and adopts the Kirchoff-Love 
hypothesis for plates by assuming that: 1) a normal to the middle surface of the 
laminate remains straight and normal which is equivalent to ignoring transverse 
shearing strains, y xz = y z = 0 and 2) that normals have constant length, e. = 0. 
The remaining strains, as derived by Jones 3 \ can be expressed 




or, equivalently, {e} = {t°} + z {*} (3.7) 

where } are the middle-surface strains and {«-} are the middle-surface 
curvatures. When equation (3.7) is substituted in equation (3. 4), the stresses in 
the &th layer can be expressed in terms of the laminate middle-surface strains 
and curvatures as 

\g x Qu Qu Q, i, £ * k x Q u Q n Qui &x 

' * ~ Q\2 Qn Qlb 'Ey ’ + Z< K y ' — Q {1 Qn Qlb ' Gy ’ AT 

T n Ql6 Q26 Qbb y" v K xy Q\b Qlb Qbb a xy 

or, equivalently, M* = |gj {e- }+ {*-} z - {«}* AT (3.8) 

where AT is the change in temperature from the stress free temperature, 7sft, to 
the operature temperature, T, which can be expressed, 

AT= T- 7sft (3.9) 

Within the scope of this dissertation, only flat, symmetric composite 
laminates with a uniform temperature are considered, and all applied forces are 
in-plane as illustrated in Figure 3.2. The geometry through the laminate 
thickness is displayed in Figure 3.3. The resultant forces, {TV}, acting on a 
laminate are obtained by integrating the stresses through the laminate thickness, 
t, and are defined as 



32 


{/V} = 


N x 

N y S- = 

N, 




o', 

<x. Wz 


(3.10) 


Equation (3.8) can be substituted in equation (3.10), and assuming the lamina 
stiffnesses to be constant within a lamina, in addition to a uniform temperature 
throughout the laminate, the relationship can be expressed as 

{/V} = [AM + -W 

where 


A=1 

A =1 

{ny 


(3.11) 

(3.12) 

(3.13) 

(3.14) 


and [A] is referred to as the laminate extensional stiffness, [B] is the laminate 
bending-extension coupling stiffness, and {N} T are the thermal forces. For 
symmetric laminates, [B] = 0, and inserting equation (3.14), equation (3.11) 









k 
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A n A l2 A l6 


£° 

X 

N 
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> = 

A\2 A 22 A 16 


£° 

-I 

6l2 022 026 

< 

a f ' 
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A | 6 A 26 A 66 J 


T° 

L *y J 

k~\ 

Q\6 Qlb Q(>6 


/V 


(z* -z k . { )AT 


or, equivalently, {N} - ( z * z *-/ (3-15) 

A =1 

Equation (3.15) can be solved for the middle-surface strains, which can be 
expressed 

{e°}= W{W}+ (316) 

where 

(317) 

A=l 

and {*} is referred to as the average laminate coefficient of thermal expansion. 
Then for in-plane loading only and [5] = 0, the middle-surface curvatures are 
{*}= o’, so substituting equation (3.16) into (3.8), the lamina stresses are solved, 
using CLT, to be 

W=0l4'{iv}t0K.Wk < 318 > 

As can be observed from equation (3.18), CLT predicts in-plane lamina stresses, 
{a}* , due to applied external in-plane loading, {/V}> and/or due to a change in 
temperature from the stress-free temperature, AT. 
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3.2 Interlaminar Stresses 

The existence of interlaminar stresses in composite laminates was first 
studied by Pipes and Pagano where they used a quasi-three-dimensional 
elasticity approach, along with finite difference numerical approximations to 
predict the interlaminar stresses, y xz , y^, and a z . 32 The problem studied is 
illustrated in Figure 3.4. For a symmetric laminate loaded axially along the ends 
x=constant, significant interlaminar stresses, y^, y^, and a z , were predicted at the 
laminate free edges, y=±b, for angle-ply laminates. At a distance of one 
laminate thickness away from the free-edge, the interlaminar stresses were 
shown to go to zero, and the stresses were equivalent to those predicted with 
classical lamination theory. Their observations included: 1) that there are no 
interlaminar stresses in unidirectional laminates, 2) interlaminar stresses are very 
high, perhaps singular, at the free edge of angle-ply laminates and 3) the 
magnitude of o z can be changed significantly by stacking sequence. 
Simplifications included: a) [ ±9] laminates exhibit only shear coupling so y« is 
the only non-zero interlaminar stress and b) [0/90] laminates exhibit only a 
Poisson mismatch between layers so yyz and a z are only non-zero interlaminar 
stresses. 

Most importantly, high positive interlaminar normal stresses may lead to 
delaminations and hence failure of the composite laminate, and high interlaminar 
shear stresses may also lead to failure and thus cannot be ignored. 
Consequently, a three-dimensional analysis is necessary to determine the 
stresses that may lead to failure in a composite laminate. 

Inclusion of resin-rich layers between plies, through the thickness of the 
laminate, was reported by Wu to affect significantly the interlaminar stresses. 3 3 
A micrograph showing a 500 times magnification through the thickness of one of 
the [0,90] s IM7/PETI5 laminate specimens is displayed in Figure 3.5. The 
micrograph shows an entire 90 degree layer located between two 0 degree plies. 
As can be observed in the micrograph, indeed there are resin-rich layers 
between lamina in the current IM7/PETI5 laminates, which should be 
incorporated into analysis models for predicting interlaminar stresses. 

3.3 Residual Thermal Stresses in Composite Laminates 

Residual thermal stresses in composite laminates arise due to curing 
orthotropic lamina in bonded layers at various orientations, from an elevated 
temperature down to the operating temperature. Residual stresses in laminates 
have been studied by a myriad of authors and is currently an active area of 
research. One of the earliest studies was conducted by Chamis in 1971, where 
he used classical lamination theoiy to predict significant in-plane residual thermal 
stresses for various laminates. 3 Temperature dependent material properties 
were included in the nonlinear analysis of Pagano and Hahn in 1977. 35 
Furthermore, a three-dimensional finite element analysis was performed by 
Chapin and Joshi in 1991. Their investigation revealed significant interlaminar 
residual thermal stresses. 
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Since the magnitudes of both in-plane and interlaminar stresses developed 
in a composite laminate due to curing have been shown to be significant, the 
extent to which they vary with a change in operating temperature is of key 
interest in this study. The change in residual stresses with a change in 
temperature has the potential of playing an important role in the change of 
strength with a change in temperature and will consequently be included in the 
analysis models developed herein. 

3.4 Validation of 3-D Finite Element Analysis Method 

A finite element analysis model is developed to replicate the analysis 
results of a previously studied problem in the literature. The purpose of this 
model is to validate the use of the current analysis tools and procedures. The 
problem of interest is a [+45,-45] s graphite epoxy laminate subject to uniform 
axial extension of 0.00005 inches, on the ends x=constant, as previously 
illustrated in Figure 3.4. Each ply is modeled as an elastic orthotropic material 
with properties: 

= 20. x 10 6 psi 
E 22 = E 33 = 2. 1 x 1 0 6 psi 
G12 = G23 = G31 = 0.85 x 1 0 6 psi 
V12 = V13 = v 2 3 = 0.21 
an = -0.5 x 10' 6 in/in°F 
0C22 = 0C33 = 23. x 10' 6 in/in°F 

The problem was studied by Pipes and Pagano in their investigation of 
interlaminar stresses. 32 They used a finite difference method along with an 
elasticity approach to solve for stresses in the y-z plane, where all stresses are 
independent of x. Their results showed extremely large interlaminar stresses, a z 
and Oxz, at the intersection of the free-edge, y=b, and at interfaces between 
layers, but these stresses approached zero within one laminate thickness from 
the free edge. The possibility of stress singularities at the intersection was 
postulated. The existence of stress singularities was further investigated by Raju 
and Crews in 1980 37 and by Raju, Whitcomb, and Goree in 1980. 8 They 
analyzed the same problem using the finite element method with 8-node, two- 
dimensional quadrilateral elements. A comparison was made between solutions 
obtained from finite difference, perturbation and finite element methods. The 
interlaminar normal stress, o z , was shown to disagree in both magnitude and 
sign. The magnitude difference was expected due to the singularity and 
approximation differences with the numerical methods; however, the sign of the 
finite element solution was concluded to be the correct sign because the other 
methods made inaccurate assumptions for the boundary conditions. 

Convergence studies conducted by Raju et al. showed the finite element 
solution to converge everywhere except very near the interface corner. 38 By 
looking as similar exact solution problems with singularities, finite element 
solutions were found to behave similarly near the singularity. Thus, the finite 
element method was concluded to be accurate except in a region involving the 
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two elements closest to the singularity, and that this region can be made 
arbitrarily small by refining the finite element model. Another study conducted by 
Reedy in 1988 showed that the singularity dominated an extremely small region 
near the free edge, and that this singularity is weak in that stresses changed by 
only 2% when the element size at the free edge was decreased by 50%. 39 
Reedy also showed that rapidly varying interlaminar stresses exist even when 
there is no singularity present, and hence that there is a real interlaminar stress 
concentration as predicted by the finite element method. The two studies by 
Raju and Reedy give validity to the use of the finite element method for predicting 
interlaminar stress concentrations very close to the free edge. As a result, in the 
present study, interlaminar stresses will be evaluated everywhere except in the 
two elements closest to the free edge interface or apparent singularity, and these 
elements will be made as small as possible, so as to improve the approximate 
values for the interlaminar stresses. 

In the current study, all finite element models are created within PATRAN, 
a commercially available finite element model generation code and then 
NASTRAN, a commercially available finite element analysis code, is used to 
perform the structural analyses. 3 10 As opposed to the quasi-three-dimensional 
models used by Raju et al. to predict stresses in the y-z plane, an actual three- 
dimensional finite element model using 20-node hexahedral elements is created 
here to solve the problem illustrated in Figure 3.4. The three-dimensional model 
is necessary for further investigation of the effect of thermal cure loading and pin- 
bearing behavior described in Chapter 5. 

Figure 3.6 shows the discretization used in the y-z plane of the first 
quadrant of Mesh 1. Mesh 2, a more refined mesh, differed in that the three 
elements closest to the free edge in Mesh 1 were divided in half in the y-z plane. 
A model containing all four in-plane quadrants of the laminate were necessary for 
a three-dimensional model with accurate boundary conditions and the lack of in- 
plane symmetry for the [+45,-45] s lay-up. Mesh 1 was chosen to closely replicate 
the refined mesh used by Raju (Figure 9(c) in Ref.3.8). Unfortunately, the exact 
nodal locations were not given by Raju. With five elements used here in the x- 
direction and 576 elements in the y-z plane of the first quadrant, there is a total 
of 1 1 ,520 hexahedral elements in the model of all four quadrants. 

The analysis results for the interlaminar stresses, o z and Oxz, are plotted in 
Figures 3.7 and 3.8, respectively. The stresses are plotted along the interface, 
z=h, for the plane x=0. These stresses replicate the interlaminar stress 
distributions reported by Raju et al. 3 8 As expected due to the singularity, Mesh 2 
predicted higher stresses in the two elements closest to the free edge than did 
Mesh 1 . Most importantly, two elements away from the free-edge of Mesh 1 , the 
magnitude of the stresses for Mesh 1 , Mesh 2, and Raju’s model all converge at 
the same value. The convergence of the interlaminar stress solution validates 
the analysis model developed here for evaluation of interlaminar stresses two 
elements from the free edge. 

Analysis results, using Mesh 1, for the in-plane stresses, o x , o y , and x xy , 
are plotted in Figure 3.9. Also shown on the figure is the classical lamination 
theory (CLT) solution for the in-plane stresses. As previously observed by other 
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authors, the analysis results for in-plane stresses can be seen to approach the 
CLT solution at a distance of one laminate thickness, 4h, from the free edge. 

The current finite element model was further utilized to investigate the 
thermal residual stresses that would be present in the laminate due to curing. A 
stress free temperature of T re f = 310°F was assumed, and the stresses were 
evaluated at an operating temperature of 70°F. The previously presented 
constant material properties were assumed, for simplicity, allowing for a linear 
static analysis. 

The in-plane stresses, a x , a y , and txy, obtained using Mesh 1 and the 
thermal curing load of A7^ -240°F are plotted in Figure 3.10. Also shown is the 
classical lamination theory (CLT) solution for the in-plane stresses. As with the 
case of axial extension loading, the analysis results for in-plane stresses also 
approach the CLT solution at a distance of approximately one laminate thickness 
from the free edge from the free edge. The figure also reveals large in-plane 
stresses near the free edges due to curing. The large o x and o y stresses appear 
to be confined to the two elements closest to the free edge and, due to the 
singularity, are consequently not valid stresses. 

The interlaminar normal, a z , and shear, a^, stresses are displayed in 
Figures 3.11 and 3.12, respectively. Similar to the case of axial extension 
previously analyzed, a stress concentration and apparent singularity are present 
at the intersection of the free-edge and interface. However, in contrast to the 
case of axial extension, the interlaminar stresses o z and a xz are positive(tensile) 
in magnitude and more severe for curing. Since the thermal residual stresses 
are also much greater in magnitude than the interlaminar stresses due to 
extensional loading alone, their significance on the changes in the state of stress 
due to a change in temperature will be further investigated within Chapters 3 and 
5. 

3.5 Three-Dimensional Laminate Analysis 

A three-dimensional analysis model of an IM7/PETI5 composite laminate 
is generated here for the purpose of establishing the laminate attributes that 
significantly affect the state of stress in the laminate. More specifically, the focus 
is on the attributes that influence the changes in stress with a change in 
temperature, and hence should be incorporated in the pin-bearing analysis 
model. The specific laminate attributes under investigation include: 1) the resin 
rich layer located between lamina in the laminate, 2) the residual thermal 
stresses that develop in a laminate due to the cure cycle, 3) the temperature 
dependent stiffness and thermal expansion coefficients for the IM7/PETI5 lamina, 
and 4) nonlinear stress-strain behavior of the IM7/PETI5 lamina. 

Due to in-plane symmetry in the [90,0]8s laminates, and the formidable 
computational size of the full three-dimensional problem, only this lay-up was 
chosen for investigation, where taking advantage of symmetry, only one-eight of 
the laminate is modeled. Two finite element geometry models of the [90,0]ss 
laminate were created, one without a resin layer and one including the resin 
layer. The finite element model solid geometry including the resin-rich layer is 
shown in Figure 3.13. The model without the resin layer differed in that there 
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were no resin layers and a ply thickness of h=0.0055 inch is used, resulting in the 
same laminate half thickness of 0.088 inch. As approximated from the 
micrograph previously shown in Figure 3.5, for the model including the resin 
layer, the ply layer thickness, h, is 0.005 inch, where the resin layer thickness is 
0.1 h. All ply layers have a thickness of h and resin layers have a thickness of 
0.1 h except the top ply has a thickness of h+0.05h and the one resin solid 
adjacent to the z=0 symmetry plane has a thickness of 0.05h. Using a value of 
h=0.005 inch, the total thickness modeled is 0.088 inch, which is half of the total 
laminate thickness of t=0.176 inch. The three symmetry planes are z=0, the 
through the thickness symmetry, and x=0 and y=0, the in-plane symmetry 
planes. 

The finite element models for Mesh 1 and the more refined mesh, Mesh 2, 
consisted of hexahedral solid elements with rectangular faces. A view of Mesh 1 
in the x-y plane is displayed in Figure 3.14, where there are 20 uniformly spaced 
elements located in the x-direction and 30 elements in the y-direction. In the y- 
direction, 19 uniformly spaced elements are between y=0 and y=0.95W, and the 
remaining eleven elements were graded with the maximum concentration located 
at the free edge, y=W. A partial view of Mesh 1 in the x-z plane, including the 
resin layer, similar to the blow-up region depicted in Figure 3.13, is displayed in 
Figure 3.15, where there are 100 element layers through the thickness of the 
model. Each ply layer had 4 elements and each resin layer had 2 elements 
except for the two resin layers between the 0° and 90° plies located second and 
third from the z=0 symmetry plane, which had four elements. A total of 38,400 
elements without the resin layer, and 60,000 elements with the resin layer, are in 
the model for Mesh 1 . Mesh 2, the more refined mesh, differed from Mesh 1 with 
35 elements in the y-direction, where the elements were made smaller than 
Mesh 1 between y=0.998W and the free edge. Also, there are 112 elements in 
the z-direction for Mesh 2 including the resin layer, which differed from Mesh 1 
with 8 elements in the first 0° ply located closest to the z=0 symmetry plane, 8 
elements in the resin layer between the first 0° and 90° ply and 4 elements per 
resin layer in the first resin layer adjacent to the z=0 symmetry plane and the 
resin layer between the first 90° and second 0° ply. A total of 48,960 elements 
without the resin layer, and 78,400 elements including the resin layer, are in the 
model for Mesh 2. 

3.5.1 Resin Laver 

To investigate the resin layer effect, the room temperature lamina 
engineering constants, as presented in Table 2.1 and Equations 2. 1-2.7, are 
used as the properties in a linear static analysis of the composite laminate. For 
the [90, 0]a s laminates, the average room temperature compressive strength of 
126 ksi, as presented in Table 2.4, is used, along with the lamina moduli, £ 11c 
and E 2 2 c , to compute an equivalent applied displacement loading of d=-0.0028 
inch for the analysis. The uniform displacement was applied on the surfaces at 
the end x=L on the finite element model as depicted in Figure 3.14. 

Analysis results for the in-plane stresses, for the finite element model 
without the resin layer, subject to the applied displacement of d = -0.0028 inch, 
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are presented in Figure 3.16, along with the Classical Lamination Theory 
solution. The most severe in-plane stress, 0 X , occurred in the 0 degree plies. 
The plots show the stresses in the 0 degree ply at the intersection of the z-0 and 
x=o symmetry planes and along the y-direction from y = 0.2 inch to the free edge 
y=0.25 inch. The in-plane stress results closely replicate the Classical 
Lamination Theory solution except in a small region near the free edge. Analysis 
results for the interlaminar stresses o z and o yz are displayed in Figure 3.17 and 
3.18, respectively. A maximum interlaminar normal stress, c z , is predicted in the 
laminate with a magnitude of approximately 7100 psi two elements from the free 
edge of Mesh 2, in the 0 degree ply along the z=0 symmetry plane. Interestingly, 
the maximum interlaminar shear stress did not occur directly at the free edge, but 
actually a very small distance away from the free edge as shown in Figure 3.18. 
A maximum of -3200 psi occurred in Mesh 2 at the interface between the 90 
degree and 0 degree ply adjacent to the z=0 symmetry plane. By comparison 
with Mesh 2 results, Mesh 1 appears to be adequate for obtaining a converged 

solution. „ . . , , 

Analysis results of in-plane stresses from the finite element model 

including the resin layers and subject to the applied displacement, were basically 
identical to the results obtained in the absence of the resin layers, however the 
interlaminar stresses are different. The interlaminar stresses predicted in the 
model including the resin layers are displayed in Figures 3.19 and 3.20^ A 
maximum interlaminar normal stress of a z = 8600 psi occurred two elements from 
the free edge of Mesh 2 along the z=0 symmetry plane. The maximum normal 
stress location was the same as obtained with the model without resin layers, 
however the magnitude was over 20% greater when the resin layers were 
included in the model. The maximum interlaminar shear stress, Cyz = -6400 psi, 
two elements from free edge in Mesh 2, occurred at the interface between the 
resin layer and 90 degree ply closest to the z=0 symmetry plane. This value was 
twice as much as the maximum interlaminar stress predicted from the model 
without the resin layers. Since there is a great difference in interlaminar stresses 
between the models with and without the resin layers, only the more detailed 
model including the resin layers will be utilized in the forthcoming analyses. 

3.5.2 Thermal Residual Cure Stresses 

The residual cure stresses in the composite laminate are investigated here 
using a linear analysis. The properties at the relevant operating temperature as 
given in Table 2.1 are used. The interlaminar stresses were the most severe 
stresses in the laminate due to curing at all three operating temperatures of 
interest -200°F, 70°F, and 350°F. Analysis results of the interlaminar stresses 
predicted at 70°F are plotted in Figure 3.21 and Figure 3.22. A maximum 
positive interlaminar normal stress of 7000 psi was predicted two elements from 
the free edge in the 90° ply at the interface of the first resin layer and 90 degree 
ply closest to z=0. A maximum interlaminar shear stress of 11500 psi was 
predicted two elements from the free edge in the resin layer at the interface of 
the first resin layer and 90 degree ply closest to z=0. Although the location of the 
maximum interlaminar stresses did not change at the different operating 
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temperatures under consideration, the magnitudes predicted did and are 
presented in Table 3.1. 


3.5.3 Temperature-Dependent Material Properties 

A nonlinear analysis of the composite laminate subject to the cure 
temperature changes is performed including temperature dependent material 
properties. A piecewise linear approximation of the material properties presented 
in Table 2.1 between the four temperatures given is incorporated in the finite 
element analysis. The temperature change loading was divided into ten 
increments with a convergence check made with a twenty increment division. 
The use of temperature dependent properties did not result in any change in the 
location of the maximum occurring interlaminar stresses and the difference in the 
magnitude in comparison to the analysis assuming constant properties was 
extremely insignificant at all operating temperatures. The maximum interlaminar 
stress results from the nonlinear analysis are presented in Table 3.1, in 
comparison to the linear analysis results. 

3.5.4 Nonlinear Material Stress-Strain Behavior 

The nonlinear elastic material stress-strain behavior was not included in 
the finite element analysis due to the following reasons, 

1) NASTRAN does not have the capability to model nonlinear elastic material 
properties for orthotropic materials, 

2) The in-plane shear stresses, X12, and interlaminar shear stress, X31, in the 
[90,0 ] 8s composite laminate being analyzed were virtually zero, and since 
only the in-plane shear modulus, G 12 , and interlaminar shear modulus, 
<J3i=Gi2, had any significant nonlinear characteristics as described in 
Section 2.3.2, the nonlinear elastic properties should not significantly alter 
the predicted stresses in the laminate analysis. 

3.5. 5 Combined Loading 

The analysis of residual cure stresses occurs from the stress free 
temperature of 460°F to the operating temperature. Then, at the operating 
temperature of interest, the applied displacement loading occurs. The case of 
cure and the case of applied displacement loading at the operating temperature 
are analyzed independently where the states of stress are subsequently 
combined. The interlaminar stresses through the thickness of the laminate, two 
elements away from the free edge are presented in Figures 3.23, 3.24, and 3.25 
at the three operating temperatures of 350°F, 70°F, and -200°F, respectively. 
Shown in each figure in the upper two plots are the ct z and Oyz stresses obtained 
independently for the case of cure loading (dashed line) and applied 
displacement loading (solid line). The bottom two plots then show the combined 
solution of cure and applied displacement loading for a z and c^. As can be 
observed from the figures at all temperatures, (except for a z at 350°F), the 
applied displacement load decreases the interlaminar residual stresses that are 



40 


present in the laminate due to cure alone. As for the c z at 350°F, the o z due to 
cure is smaller than o z due to the applied displacement, and consequently the 
net o z in the combined case is larger than o z from cure alone. However, the 
combined maximum magnitude is only 2700 psi, which is much less the o z 
magnitudes of approximately 7000 psi observed due to cure only at the 70 F and 

350°F operating temperatures. . _ 

The only significant in-plane stress in the analysis model is o x in the u 

degree plies, where i xy was virtually zero and o y was much less than the y- 
direction strength, F 22 c, in the 90 degree plies. Table 3.2 shows o x due to cure 
only, applied displacement only, and the combined case of cure and applied 
displacement loading. 

3.6 Discussion and Conclusions 

From the analysis results of the IM7/PETI5 [90,0] 8s composite laminate 
presented in Section 3.5, the significant laminate attributes that need to be 
included in the analysis of the laminate include: 1) the residual stresses due to 
cure and 2) the resin layer located between the lamina in the laminate. 
Moreover, analyses for this particular laminate at the temperatures being 
considered, establish that the temperature dependent material properties do not 
require a nonlinear analysis to provide an accurate stress solution. Just 
modeling the laminate assuming the properties to be constant at the operating 
temperature of interest as presented in Table 2.1 is sufficient for predicting 
stresses accurately. Also, due to the in-plane shear and interlaminar shear, T13, 
being virtually zero in the laminate under consideration, the nonlinear elastic 
characterization of the material is assumed not to be necessary to analyze the 

laminate behavior. . , . 

From the observation on the effect of an applied displacement load 

decreasing the interlaminar stresses present in the laminate due to cure alone, it 
is concluded that failure cannot be attributed to the interlaminar stresses in the 
composite laminate being studied when a compressive load is applied. The 
strength changes of F 11c and Fx* for the [90,0] 8s laminate, from Tables 2.1 and 
2.4, respectively, are presented in Table 3.3 along with the change in the o x 
residual cure stress in the 0 degree plies. Unfortunately, due to the small 
changes in residual cure stresses, less than 15 ksi, and a scatter of up to 15 ksi 
in the average values of F 1t0 and F xxc data presented in Table 2.1, the effect of 
residual cure stresses and changes in lamina compared to laminate strength 
cannot be clearly perceived. One can conclude that the increase in compressive 
strength of the [90,0] 8s laminate due to a decrease in temperature is most likely 
governed by the increase is lamina strength, Fn c , with a decrease in temperature 
in addition to possible effects with the changes in residual stress lessening the 
strength change for the laminate as compared to the lamina. Although 
interlaminar stresses do not appear to cause failure in a compression specimen, 
investigating the three-dimensional state of stress including thermal residual 
stresses was necessary to make the determination. Similarly, the three- 
dimensional state of stress in a pin-bearing specimen must be studied to 
determine the role of interlaminar stresses in the pin-bearing problem. 
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Table 3.1 Comparison of Predicted Maximum Residual Cure Stresses Two 
Elements Away from Free Edge. 


Temperature 

(°F) 

Linear Analysis 
Constant Properties 

w — 

Nonlinear Analysis 
Temperature-Dependent Properties 

O z (psi) 

Oy z (pSi) 

Gz(psi) 

Oyz(pSi) 

350 

1521. 

2439. 

1524. 

2439. 

70 

7081. 

11467. 

7174. 

11491. 

-200 

7691. 

12375. 

7782. 

12408. 


Table 3.2 Maximum In-Plane Stress, a x , Predicted in the Laminate. 


Temperature 

(°F) 

o x (ksi) 

Cure Only 

Applied 

Displacement 

Combined 

350 

-3. 

-230. 

-233. 

70 

-15. 

-231. 

-246. 

-200 

-16. 

-231. 

-247. 


Table 3.3 Comparison of Changes in Lamina Strength, Fn c , [90,0]8 S Laminate 
Strength, Fxxc, and o x Cure Stress in 0 degree Ply with a Change in Temperature. 


W Vi . j-j ... • , ■ AAU — ' ■ 

Temperature 

Change 

AFn c Lamina 
(ksi) 

AF XXC Laminate 
(ksi) 

Cure Stress, Aa x 
(ksi) 

350°F to 70°F 

40 

31 

-12 

350°F to -200°F 

53 

49 

-13 
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Figure 3.1 . Illustration of transformation between the material and laminate 
coordinate systems. 



Figure 3.2. Composite laminate subject to in-plane applied forces. 



Figure 3.3. Geometry through the thickness of an N-layered laminate. 
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w xy 

c) 3-D stress components. 


Figure 3.4. Illustration of the problem of a four-ply laminate under uniform axial 
extension and the resultant stress components, Ref.3.2. 
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Figure 3.5. Micrograph showing a 500 times magnification of a portion through 
the thickness of a (90,0) s IM7/PETI5 laminate. 
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Figure 3.6. Finite element model discretization of Mesh 1 in first quadrant of the 
y-z plane. 
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Figure 3.7. The interlaminar normal stress, o Zt along the interface, z=h. 
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Figure 3.8. The interlaminar shear stress, Oxz. along the interface, z-h. 
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Figure 3.12. 



Figure 3.13. Finite 
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Figure 3.14. Finite element discretization in the x-y plane for Mesh 1 . 
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Figure 3.16. In-plane stresses in 0° ply of [90,0] 8s laminate, subject to 
applied displacement, plotted along the z=0 and x=0 symmetry planes. 
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Figure 3.20. Interlaminar shear stress in [90,0] 8s laminate including resin 
layers, subject to applied displacement, plotted along x=0 and the interface of the 
first resin layer and 90 degree ply closest to z=0. 
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Figure 3.22. Interlaminar shear stress in [90,0] 8s laminate at 70°F due to curing 
only, plotted along x = 0 in the resin layer at the interface of the first resin layer 
and 90 degree ply closest to z=0. 
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Displacement Cure 



o z , psi G yz , psi 


combined combined 



o z , psi psi 

Figure 3.23. Interlaminar stresses through the thickness of the laminate two 
elements away from free edge at T=350°F. 
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Figure 3.24. Interlaminar stresses through the thickness of the laminate two 
elements away from free edge at T=70°F. 
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Figure 3.25. Interlaminar stresses through the thickness of the laminate two 
elements away from free edge at T=-200°F. 
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Chapter 4 
Pin-Bearing Tests 


An experimental investigation on the effect of temperature changes on the 
pin-bearing strength is pursued to establish the extent of bearing strength 
changes with a change in temperature for the IM7/PETI5 material system. 
Several lay-ups were investigated to observe the effect of lay-up on the extent of 
the strength change and to establish a possible optimal lay-up for maximum 
strength at varying temperatures. The test data for the [90,0] 8s lay-up is utilized 
for validation of the computational model through qualitative comparison of 
experimental strength changes with temperature compared to changes in the 
state of stress at the bolt hole of the pin-bearing computational model with a 
change in temperature. Within this chapter the test procedure is described, 
followed by test configuration and mechanical operation, the test specimen 
details, and finally the test results. 

4.1 Test Procedure 

The pin-bearing tests were conducted within the guidelines of the 
American Society of Testing and Materials (ASTM) standard for determining the 
bearing response of a polymer matrix composite laminate. 41 The ASTM test 
method was also consistent with the recommendations of MIL-HDBK-17. The 
double shear tensile loading procedure was chosen to provide pin bearing loads 
on a flat, rectangular composite laminate coupon with a centerline hole. The 
thickness of the composite at the hole boundary and the hole diameter are 
measured prior to testing. The composite laminate coupon is loaded at the hole 
in bearing as illustrated in Figure 1.4. The bearing load is applied through a 
close tolerance pin that is reacted in double shear by the fixture, creating a 
symmetric loading condition through the thickness of the composite laminate 
hole. The load is applied to the composite laminate specimen while both the 
load and hole deformation are monitored. The test is terminated when a 
maximum load has clearly been reached, and the load is removed immediately 
so as not to mask the initiation failure characteristics with additional damage. 
The bearing stress versus bearing strain is plotted and the failure mode is noted. 
The ultimate bearing strength is determined from the maximum load and the 
ultimate bearing strain is determined after correcting the strain data for the new 
effective origin as will be described later is Section 4.3. 

4.2 Test Configuration and Operation 

A schematic drawing of the test configuration is shown in Figure 4.1. A 
similar configuration was used by Crews in 1986 and later modified by Wright in 
1997. 4344 The test components shown in the figure include the test specimen, 
the pin for applying the bearing load to the test specimen, the extensometer and 
measurement wire for measuring the bolt hole elongation, and the clevis and grip 
which are the fixtures connecting the specimen to the test machine. The pin- 
bearing load is reacted in double shear by the clevis which is attached to the 
load frame while the end of the specimen constrained by the bolted grip is held 
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fixed. Photographs of the test configuration in the 22kip test machine, 
surrounded by the environmental test chamber, are shown in Figure 4.2. 

The pin used for applying the bearing load is made of Vascomax C-300 
steel. 44 The pin material was chosen for its high strength characteristics and 
ability to transfer high bearing loads without the possibility of pin failure during 
testing. The manufacturer’s available pin properties of interest included, 

d = 0.249 in 

E p = 27.5 x 10 6 psi 

oc p = 5.6 x 10' 6 in/in/F 

F cy = 317. ksi 

F-rv = 290. ksi at T=70°F and F^ = 245. ksi at T = 600°F 

The pin that goes through the specimen bolt hole is also placed through 
clearance fit bolt holes on both sides of a steel clevis in a double shear 
arrangement. An attachment rod connects the clevis with the load frame. The 
lower end of the specimen is held constrained with a steel bolted grip which is 
also connected to the test machine with an attachment rod as can be seen in the 
photograph in Figure 4.2. A 30 ft-lb torque was used for all six grip bolts. 

The extensometer and measurement wire are used to measure the bolt 
hole elongation during testing. The measurement wire is made of 0.060 inch 
diameter stainless steel wire. A high temperature epoxy adhesive is used to 
fastener one end of the wire in the slot machined to fit the wire in the unloaded, 
bottom portion of the specimen bolt hole, and the other end of the wire is 
attached to the extensometer. The extensometer measures the relative 
displacement between its two arms. One arm is attached to the end of the clevis 
as shown in Figure 4.1 where the other arm attaches to the measurement wire. 
Assuming negligible pin bending, the displacement of the loaded top of the bolt 
hole is measured relative to the unloaded end where the measurement wire is 
fastened, giving the bolt hole elongation. 

The environmental test chamber with the door open can be seen in the 
Figure 4.2 photograph, wherein the specimen is either convectively heated or 
cooled to the desired test temperature prior to applying the mechanical bearing 
load. For the -200°F and 350°F temperature tests, three thermocouples are 
monitored, one on each side of the specimen nearly adjacent to the bolt hole and 
one located on the clevis. After approximately one hour, when all thermocouple 
temperatures remain steady within two degrees of the desired operating 
temperature, the mechanical loading is applied to the specimen. The test 
machine is operated by displacing the clevis at a rate of 0.004 in/min while the 
mechanical load and bolt hole elongation are measured and recorded each 
second. The test is terminated at the onset of failure as observed by an obvious 
drop in the load. 

4.3 Test Specimens 

The composite laminate test specimens for the pin-bearing tests were 
designed with 32 plies to transfer large loads and hence be representative of a 
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joint in a high speed aircraft wing or launch vehicle intertank. The objective was 
also to design the specimens with sufficient edge, e, and width, w, distances to 
fail in the bearing mode for all test temperatures and eliminate premature failure 
by shearout or tension, seen in Figure 1 .6. The test specimen geometry is shown 
in Figure 4.3. The length, L= 8 inches, was chosen to provide adequate gripping 
area to prevent any slipping on the fixed end of the specimen. The bolt holes 
were drilled in the specimens by first drilling under the specified diameter, D, and 
then reaming to size. A small hole of diameter 0.063 inch was also drilled 
adjacent to the bolt hole as shown, where the material between the small hole 
and bolt hole was then filed away to allow insertion of the measurement wire 
through the bolt hole. A specimen width, W = 2 inches, was used for all 
specimens except the specimens fabricated from the [90] 32 lay-up, where W- 4 
inches, in an attempt to avoid a tension failure in these specimens with low 
tensile strength. A photograph of a 2. inch wide pin bearing test specimen is 
provided in Figure 4.4. 

Pin-bearing specimens were made using five different laminate lay-ups as 
defined in Table 4.1. For each lay-up, three 12 in by 12 inch panels were 
fabricated where five replicate specimens were made from each panel. A total of 
105 pin-bearing test specimens were fabricated for the pin bearing experiments. 


4.4 Test Results 

The pin-bearing experiments were conducted to determine the pin-bearing 
strength for all lay-ups tested at the three operating temperatures of 350°F, 70°F, 
and -200°F. The pin-bearing stress - bearing strain curve for one of the [90,0] 8s 
laminates tested at 70°F is shown in Figure 4.5. This curve was typical for all 
specimens tested successfully. The pin-bearing stress is calculated as 



(4.1) 


where P is the measured load transmitted , D is the bolt hole diameter, and t is 
the bolt hole thickness. The pin-bearing strain is calculated as 


E 


br 


b_ 

D 


(4.2) 


where 8 is the bolt hole elongation measured by the extensometer. The ultimate 
bearing strength is then determined as 


p 

p max 

bru Dt 


where P mm is the maximum load attained during testing, 
strain is determined as 

^bru ^brrmax ^0 


(4.3) 

The ultimate bearing 

(4.4) 

where e is the bearing strain at the maximum load, and e 0 is the bearing strain 
at the effective origin. In accord with the ASTM standard, an effective origin is 
computed from the initial linear portion of the stress-strain curve as illustrated in 
Figure 4.5. The initial nonlinearity is typical of a bearing test and attributable to 
initial joint translations and straightening. 4 1 The effective origin is computed from 
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the chord stiffness line intersection with the bearing strain axis. Similar to e bru , the 
bearing strain, e^, is determined from the bearing yield strength, F bry . 

For the [90,0] Ss specimens tested, in addition to the ultimate bearing 
strength, the ultimate bearing strain, bearing chord modulus, and bearing yield 
strength and bearing yield strain were also computed for comparison with the 
analysis models developed in Chapter 5. The chord modulus, E br , is calculated 
from the slope of the linear portion of the stress-strain curve, where the curve 
between c br = 25 ksi and o br = 40 ksi is used in this report. For several 
specimens, the extensometer did not function properly during testing, but the 
ultimate bearing stress was still determinable. The test data for all specimens of 
Layup 1, the [90, 0]^ laminate, is presented in Table 4.2. Also calculated are the 
average values for F bnj , , and £ br . The ultimate bearing strength, F bru , data for 
Layup 2, 3, 4, 5, 0, and 90 is presented in Tables 4.3- 4.8, respectively. Also 
presented in Tables 4. 3-4. 8 are the failure modes for each specimen. 

Excluding the [0] 32 and [90] 32 lay-ups, all specimens tested failed in the 
bearing mode as illustrated in Figure 1 .6, which was the desired failure mode for 
these tests. Some of the [0] 32 specimens failed in the shearout mode, and the 
number at -200°F that failed in shearout was greater than at 70°F, but no 
specimens tested at 350°F failed in shearout. This indicates that decreasing 
temperature increases the likelihood of the shearout mode, whereas as 
temperature increases, bearing strength decreases and the bearing failure mode 
becomes more likely. Although the [90] 32 specimens were made twice as wide as 
all other specimen, that was still not sufficient and most [90] 32 specimens failed 
in the tension mode, which would be expected for the [90]^ specimens with low 
tensile strength. 

The average ultimate bearing strength data for all specimens is plotted in 
Figure 4.6. Error bars show the range of values associated with the plotted 
average value. From this curve, it is obvious that for specimens that failed in the 
bearing mode, the bearing strength increases with decreasing temperature for all 
the lay-ups tested. Although the average bearing strength of Layup 3 was the 
greatest at all temperatures, due to the range in values at each temperature, it is 
not clear that this [+45,-45,0,90] 4S quasi-isotropic lay-up is necessarily the 
optimum lay-up. Layup 1, with a [90,0] 8s lay-up and Layup 4 with a [+45,0,- 
45,90] 4S lay-up also faired well with the range of bearing strengths measured 
overlapping Layup 3. For the [90] 32 lay-up which failed in the tension mode, the 
strength remained almost constant as does the tensile strength of the material 
with varying temperature. 

Further evaluation of some failed specimens included having the 
specimens carefully cut along the bearing load line for observation under a 
microscope. To observe internal damage, the specimens were cut along the 
dashed lines illustrated in Figure 4.7. Specimens where the external damage 
was minimally visible from Layup 1 , 3, 4, and 5 were chosen to observe internal 
failure initiation characteristics. For Layup 1 , the [90,0] 8s laminate, specimens c- 
1, b-5 and a-4, as identified in Table 4.2, were evaluated. All three specimens 
have very similar internal failure characteristics. A series of three micrographs of 
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specimen c-1, tested at -200°F, with a nearly 50 times magnification, is 
presented in Figure 4.8. Extreme damage is evident in the outermost plies, with 
a transverse shear crack emanating from that region. These failure 
characteristics were also observed within specimens b-5 and a-4. The 
micrographs indicate the possibility of failure initiating in the [90,0] 8s specimens 
due to the interlaminar shear stress at all test temperatures. These observations 
will be compared with pin-bearing analysis results in Chapter 5, to determine if 
there is any correlation with significant stresses predicted with the analysis 
model. 

Several observations can be made in comparing the average pin-bearing 
ultimate strengths from Tables 4.2-4.6 to the corresponding laminate strengths 
presented in Table 2.1. First, pin-bearing strength is observed to be significantly 
lower than the corresponding laminate compressive strength in most cases. 
Only for Layup 3 and 4 at -200°F pin-bearing strength was greater and for Layup 
4 at 70°F pin-bearing strength and laminate strength were equivalent. Also, the 
extent of pin-bearing strength degradation in the temperature range from 70°F to 
350°F is observed to be lower than the laminate strength degradation for Layup 1 
and 5, but the opposite is true for all other lay-ups in Table 2.1, where the pin- 
bearing strength degradation was greater than the laminate strength 
degradation. With a temperature decrease in the temperature range from 70°F 
to -200°F, the pin-bearing strength increase was greater than the laminate 
strength increase for all lay-ups. Although both laminate strength and pin- 
bearing strength similarly decrease with an increase in temperature, the extent of 
the laminate strength change differs from pin-bearing strength change. In 
general, the extent of both laminate strength change and pin-bearing strength 
change with a change in temperature is dependent on the lay-up of the 
IM7/PETI5 material and the temperature range. 
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Table 4.2 Pin-Bearin g Test Data for Layup 1 - [90,0] 8s . 

T = -200°F 

t I D I P ^ R I e I E F 

max bru °bru br 1 bry 

(in) (in) (lb) (ksi) (in/in) (Msi) (ksi) < 




0-1 

.168 

.249 

4930 

117.8 

— 

— 


c-2 

.171 

.250 

5005 

117.1 

.0475 

2.88 

62. 

c-3 

171 

.249 

5417 

127.2 

.0510 

2.97 

62. 

c-4 

.172 

.249 

5094 

118.9 

.0560 

2.51 

64. 

c-5 

.171 

.249 

5412 

127.1 

.0610 

2.50 

65. 

Average 

Value 




121.6 

.0539 

2.72 

63. 




T 

’ = 70°F 




1 

t 

(in) 

D 

(in) 

P max 

(lb) 

P bru 

(ksi) 

Ebru 

(in/in) 

E* 

(Msi) 

P bry 

(ksi) 

a-1 

.165 

.252 

3241 

77.9 

.0471 

1.91 

50. 

a-2 

.175 

.249 

3881 

88.9 

.0517 

2.14 

48. 

a-3 

.173 

.249 

3978 

91.9 

.0504 

2.14 

48. 

b-4 

.169 

.249 

3697 

87.7 

.0400 

2.68 

47. 

b-5 

.165 

.249 

3800 

92.1 

.0470 

2.63 

47. 

Average 

Value 




87.9 

.0472 

2.31 

48. 




T 

= 350°F 




1 

t 

(in) 

D 

(in) 

p ™* 

(lb) 

P bnj 

(ksi) 

^bnj 

(in/in) 

E* 

(Msi) 

P bry 

(ksi) 

b-1 

.169 

.249 

2673 

63.5 

.0466 

1.40 

41. 

b-2 

.166 

.249 

2467 

59.7 

— 

— 

— 

b-3 

.169 

.249 

2721 

64.7 

— 

— 

— 

a-4 

.172 

.249 

2623 

61.2 

.0459 

1.32 

45. 

a-5 

.168 

.249 

2646 

63.3 

.0477 

1.48 

41. 

Average 

Value 




62.5 

.0467 

1.40 

42. 


ll 






Failure Mode: All Layup 1- [90,0] Ss specimens failed in the Bearing Mode 





T able 4.3 Pin-Bearing Test Data for Layup 2 - [45,-45; 

T = -200°F 


8s* 



t 

D 

p ma * 

^bru 

Failure 

2 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

c-1 

.167 

.249 

5017 

121 

Bearing 

c-2 

.178 

.249 

— 

— 

— 

c-3 

.180 

.249 

5466 

122 

Bearing 

c-4 

.179 

.249 

5413 

121 

Bearing 

c-5 

.174 

.249 

5076 

117 

Bearing 

Average 




120 


Value 







T = 

70°F 




t 

D 

P ma* 

^bru 

Failure 

2 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

a-1 

.167 

.249 

3616 

87 

Bearing 

a-2 

.176 

.249 

4093 

93 

Bearing 

a-3 

.177 

.249 

4175 

95 

Bearing 

a-4 

.176 

.249 

3895 

89 

Bearing 

a-5 

.174 

.249 

3462 

80 

Bearing 

Average 




89 


Value 







T = 

350°F 




t 

D 

P max 


Failure 

2 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

b-1 

.170 

.249 

2499 

59 

Bearing 

b-2 

.187 

.249 

2750 

59 

Bearing 

b-3 

.187 

.249 

2487 

53 

Bearing 

b-4 

.187 

.249 

2830 

61 

Bearing 

b-5 

.179 

.249 

2221 

50 

Bearing 

Average 




56 


Value 







T abl e 4.4 Pin-Bearing Test Data for Layup 3 - [+45,-45,0,9 01,. 

T = -200°F I 


3 

t 

(in) 

D 

(in) 

p ™, 

(lb) 

P bru 

(ksi) 

Failure 

Mode 

c-1 

.177 

.249 

5450 

124 

Bearing 

c-2 

.179 

.250 

5680 

127 

Bearing 

c-3 

.177 

.250 

5408 

122 

Bearing 

c-4 

.178 

.250 

5647 

128 

Bearing 

c-5 

.179 

.249 

5534 

124 

Bearing 

Average 

Value 




125 


T = 70°F 

3 

t 

(in) 

D 

(in) 

P max 

(lb) 

P bru 

(ksi) 

Failure 

Mode 

a-1 

.166 

.249 

4136 

100 

Bearing 

a-2 

.177 

.249 

4440 

101 

Bearing 

a-3 

.178 

.249 

4146 

94 

Bearing 

a-4 

.177 

.249 

4557 

101 

Bearing 

a-5 

.181 

.249 

3573 

79 

Bearing 

Average 

Value 




95 


T = 350°F 

3 

t 

(in) 

D 

(in) 

P max 

(lb) 

(ksi) 

Failure 

Mode 

b-1 

.169 

.249 

2777 

57 

Bearing 

b-2 

.170 

.249 

2829 

67 

Bearing 

b-3 

.173 

.250 

2671 

62 

Bearing 

b-4 

.173 

.249 

2758 

64 

Bearing 

b-5 

.179 

.249 

2811 

63 

Bearing 

Average 

Value 




63 





Tabl e 4.5 Pin-Bearing Test Data for Layup 4 - [+45,0,-45,9 0]^ 

T = -200°F 



t 

D 

P ma, 

^"bru 

Failure 

4 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

c-1 

.174 

.250 

5394 

124 

Bearing 

c-2 

.184 

.249 

5854 

128 

Bearing 

c-3 

.185 

.250 

5708 

123 

Bearing 

c-4 

.184 

.249 

5696 

124 

Bearing 

c-5 

.184 

.250 

5759 

125 

Bearing 

Average 




125 


Value 







T = 

70°F 




t 

D 

P max 

^~bru 

Failure 

4 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

a-1 

.164 

.249 

3685 

90 

Bearing 

a-2 

.180 

.249 

4203 

94 

Bearing 

a-3 

.182 

.249 

4300 

95 

Bearing 

a-4 

.170 

.250 

4205 

99 

Bearing 

a-5 

.158 

.250 

3531 

89 

Bearing 

Average 




93 


Value 







T = 

350°F 




t 

D 

P max 

F*. 

Failure 

4 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

b-1 

.171 

.250 

2168 

51 

Bearing 

b-2 

.172 

.250 

2613 

61 

Bearing 

b-3 

.174 

.247 

2714 

63 

Bearing 

b-4 

.181 

.249 

2324 

52 

Bearing 

b-5 

.177 

.249 

1775 

40 

Bearing 

Average 




53 


Value 








Table 4.6 Pin-Bearing Test Data for Layup 5 - 

r+45,0„-45,0„+45,0„-45,0„+45,90„-451,. 




T = 

-200T 




t 

D 


^”bru 

Failure 

5 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

c-1 

.168 

.250 

4511 

107 

Bearing 

o-2 

.180 

.249 

5141 

115 

Bearing 

c-3 

.181 

.250 

4916 

109 

Bearing 

c-4 

.181 

.251 

4976 

110 

Bearing 

c-5 

.172 

.249 

5056 

118 

Bearing 

Average 

Value 




112 




T : 

= 70°F 




t 

D 

P max 


Failure 

5 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

a-1 

.167 

.249 

3477 

84 

Bearing 

a-2 

.178 

.249 

3914 

88 

Bearing 

a-3 

.184 

.250 

3986 

87 

Bearing 

a-4 

.176 

.249 

3980 

91 

Bearing 

a-5 

.162 

.250 

3426 

85 

Bearing 

Average 




87 


Value 







T = 

350°F 




t 

D 

P max 

F[>ru 

Failure 

5 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

b-1 

.175 

.249 

2307 

53 

Bearing 

b-2 

.185 

.249 

2340 

51 

Bearing 

b-3 

.186 

.250 

2593 

56 

Bearing 

b-4 

.180 

.250 

2370 

53 

Bearing 

b-5 

.169 

.250 

2338 

55 

Bearing 

Average 




54 


Value 







Table 4.7 Pin-Bearing Test Data for Layup 0 - [OL 




T = 

-200°F 




t 

D 

P ma* 

^bru 

Failure 

0 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

c-1 

.156 

.250 

3780 

97 

Shearout 

c-2 ' 

.176 

.250 

4252 

97 

Shearout 

c-3 

.182 

.249 

4177 

92 

Shearout 

c-4 

.181 

.250 

4348 

96 

Shearout 

c-5 

.174 

.250 

3514 

80 

Bearing 

Average 




93 


Value 







T = 

= 70°F 




t 

D 

P ma , 


Failure 

0 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

a-1 

.177 

.250 

2964 

67 

Bearing 

a-2 

.185 

.250 

3238 

70 

Shearout 

a-3 

.182 

.250 

3320 

73 

Bearing 

a-4 

.179 

.250 

3197 

71 

Shearout 

a-5 

.162 

.249 

2762 

68 

Shearout 

Average 




70 


Value 







T = 

= 350°F 




t 

D 

P max 


Failure 

0 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

b-1 

.160 

.249 

1535 

39 

Bearing 

b-2 

.177 

.250 

1539 

35 

Bearing 

b-3 

.184 

.251 

1635 

35 

Bearing 

b-4 

.182 

.250 

1717 

37 

Bearing 

b-5 

.176 

.250 

1604 

36 

Bearing 

Average 




36 


Value 








Table 4.8 Pin-Bearing Test Data for Layup 90 - f901,. 




T = 

-200°F 




t 

D 

P ma* 

^bru 

Failure 

90 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

c-1 

.175 

.249 

969 

22 

Tension 

c-2 

.175 

.249 

1009 

23 

Tension 

c-3 

.175 

.249 

962 

22 

Tension 

c-4 

.186 

.249 

1118 

24 

Tension 

c-5 

.188 

.249 

1344 

29 

Tension 

Average 

Value 




24 




T : 

= 70°F 




t 

D 

p ™. 

bru 

Failure 

90 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

a-1 

.183 

.249 

1198 

26 

Tension 

a-2 

.184 

.249 

1190 

26 

Tension 

a-3 

.181 

.249 

1233 

27 

Bearing 

a-4 

.184 

.249 

1153 

25 

Tension 

a-5 

^ .185 

.249 

1080 

23 

Bearing 

Average 

Value 




25 




T = 

350°F 




t 

D 

Pma, 

P bru 

Failure 

90 

(in) 

(in) 

(lb) 

(ksi) 

Mode 

b-1 

.184 

.249 

1171 

26 

Bearing 

c-2 

.183 

.250 

1217 

27 

Tension 

b-3 

.181 

.249 

1207 

27 

Tension 

b-4 

.185 

.249 

1151 

25 

Tension 

b-5 

.186 

.250 

1205 

26 

Tension 

Average 

Value 




26 




Clevis 



Grip 

Figure 4.1 . Schematic drawing of pin-bearing test configuration. 




a) 


b) close-up 

Figure 4.2. Photographs of test set-up in test machine. 
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L = 8. in 

◄ - > 

e = 1.5 in 
< ► 



Figure 4.3. Pin-bearing test specimen geometry. 
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Figure 4.4. Photograph of pin-bearing test specimen. 






86 


Tension 

Plane 

r 



* 


Bearing 

Plane 


Bearing 

Surface 



Figure 4.7. Illustration of failed specimen cut for observation of internal damage. 
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Figure 4.8. Series of three micrographs showing damage through the thickness 
at the bolt hole of specimen 1-c of the [90,0] 8s specimens tested at -200°F. 
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Chapter 5 

Computational Study 

Computational models are developed to study the effect of temperature 
changes on the state of stress in a composite laminate subject to pin-bearing 
loads. The significance and role of each basic change, AS ,A M, AS, as defined by 
the hypothesis of this dissertation (see Chapter 1 .3), is investigated. The finite 
element model of the composite laminate containing the bolt hole is first 
presented. Two preliminary studies are then introduced, one aimed at 
understanding the effect of temperature change on the state of residual stress in 
the vicinity of the bolt hole, AS, and another one for understanding the effect and 
significance of temperature change on the fastener fit, AS. The pin-bearing 
problem is then studied by first applying the pin-bearing load using the material 
properties at the three different operating temperatures. A solution combining 
residual stresses with pin-bearing stress results is presented. Finally, 
comparisons are made between computationally and experimentally determined 
results. 

5.1 Finite Element Model 

Following the guidelines developed in Chapter 3, the finite element 
method is used to perform three-dimensional analyses of a composite laminate 
containing a bolt hole. The [90,0] 8s laminate is analyzed where in-plane 
symmetry as well as out-of-plane symmetry allow the computational model to be 
reduced to one-fourth of the actual configuration. This section presents model 
attributes that are consistent for all analyses to be conducted within this chapter. 
These model attributes include the geometry, finite element mesh, and material 
properties. 

Geometry. The model geometry is illustrated in Figure 5.1. Symmetry planes 
exist about z=0, the laminate thickness symmetry, and about y=0, the in-plane 
symmetry. Although the geometry is symmetric about x=0, the pin-bearing 
loading condition is not and hence the model cannot be further reduced. The in- 
plane size of the laminate modeled was 0.5 inch wide and 2. inches in length. 
Although the pin-bearing test specimens were 2. inches wide and 8. inches in 
length, to reduce the computational size, the laminate size modeled was 
reduced. Through evaluation of analysis results forthcoming, the in-plane size of 
the laminate modeled proved to be sufficiently large, so that laminate size did not 
significantly alter the state of stress in the vicinity of the bolt hole. The through 
the thickness laminate geometry includes both ply layers and resin layers, where 
a total of 32 solid layers are contained in the model, wherein 8 layers are 0 
degree plies, 8 layers are 90 degree plies, and 16 layers are resin layers. This is 
the same thickness laminate geometry modeled in Chapter 3 and depicted in 
Figure 3.13. Once again the ply thickness is h = 0.005 inch and the resin layer 
thickness is O.lh, where only half the resin layer is modeled adjacent to the z=0 
symmetry plane and the top 90 degree ply has a thickness of 0.00525 inch. The 
total laminate thickness modeled was 0.088 inch, which is half the actual 
laminate thickness of 0.176 inch. 
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Finite Elements. The finite element model was generated using 8-node 
hexahedral solid elements. A convergence study, using the room temperature 
thermal cure loading condition stress results, was conducted to determine an 
optimum mesh. The objective was to refine a mesh until all stress components 
were converged (when compared to a more refined mesh) two elements away 
from the free edge bolt hole boundary, but also so that the mesh was not more 
refined than necessary. Although less refinement may have been achievable 
with further study, an acceptable mesh was defined which contained a total of 
63,498 elements. Figure 5.2 shows the finite element mesh in the x-y plane. 
Figure 5.2(a) shows the complete finite element model and Figure 5.2(b) displays 
an enlarged view showing just a 0.25 inch by 0.25 inch section to the right of the 
bolt hole on the z=0 plane. This section is where the largest stresses are 
expected in the pin-bearing model. Figure 5.2(b) also shows the origin, a 
rectangular coordinate system and a cylindrical coordinate system at the center 
of the bolt hole. 

The two layers of elements closest to the bolt hole boundary were only 
1.67x1 O' 5 inch wide in the radial direction. The analysis results obtained in the 
two elements closest to the bolt hole boundary, within which the free edge 
singularity is contained, will not be presented. Alternatively, analysis results are 
presented at a distance of 3.3 x 10 5 inch in away from the bolt hole boundary 
where r = 0.125033 inch; hence all stresses presented are assumed to be valid 
computational stresses. 

At the bolt hole boundary, a total of 108 layers of elements were located 
through the thickness of the laminate in the z-direction. The first resin solid layer 
adjacent to the z=0 symmetry plane and the next two plies (a zero degree ply 
and 90 degree ply) were modeled with five layers of uniformly spaced elements 
each. The remaining solid layers were discretized each with three layers of 
uniformly spaced elements, except the three outermost layers closest to z=0.088 
inch. The three outermost layers were discretized each with five layers of graded 
elements so that the elements were smaller closest to the interfaces between 
layers. Between r = .1251 inch and r=.145 inch, the through the thickness mesh 
was transitioned to only one element per solid layer through the thickness of the 
laminate. 

Material Properties. The material properties used in the analyses were assumed 
to be constant at the applied operating temperature. Through comparison with 
nonlinear temperature-dependent property analyses conducted in Chapter 3, the 
constant property assumption was determined to be sufficient for accurately 
predicting stresses at the operating temperatures being studied. For the 
operating temperatures of 350°F, 70°F, and -200°F, the material properties, as 
given in Table 2.1, were used. Since compression was the primary loading 
condition, compressive stresses were assumed to dominate. Hence for 
simplicity, and although tensile stresses were predicted with the analysis model, 
only the compressive moduli, E 1ic and E 22 c, are used in the analysis models. 
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5.2 Residual Stress Concentrations at Bolt Holes 

The state of residual stress due to curing a composite laminate is an initial 
state of stress, and it is present in the laminate containing a bolt hole prior to 
applying the pin-bearing load. The study presented focuses on predicting the 
changes in the state of thermal residual stress, AS, in the vicinity of the bolt hole 
due to the change in temperature. The change in stress, AS, for the current 
problem is a function of both the change in temperature, AT, and the change in 
material properties of the constituent materials, AM. To predict the change in the 
state of residual stress accurately from room temperature to an elevated or 
reduced temperature, knowledge of the residual stress concentrations due to 
curing present in the laminate at room temperature is required. 

In this study, the composite laminates were first cured, where 
subsequently, at room temperature, the bolt hole was drilled and reamed to size 
in the composite laminate having an initial state of residual cure stresses. 
However, the finite element analysis for evaluation of residual stresses at a room 
temperature of 70°F is performed on a model containing the bolt hole at the 
stress free temperature of 460°F, where a temperature change to room 
temperature is applied as the thermal load. Moreover, the state of residual stress 
in the vicinity of the bolt hole is assumed here to be altered insignificantly by 
modeling the bolt hole in the laminate at the stress free temperature, and then 
curing to room temperature. 

In the three previous published studies on thermal residual stresses in 
composite laminates containing a bolt hole the same modeling assumption was 
used where the bolt hole was assumed present throughout the cure cycle. In the 
first study 51 , three-dimensional finite element analyses of 6-ply symmetric 
laminates were performed on a model with a very crude mesh containing only 32 
elements with only three elements through the thickness. Only the interlaminar 
normal stress result was presented at the bolt hole boundary. Accordingly, the 
interlaminar stress result presented, being at the free edge bolt hole boundary, 
must be invalid as discussed earlier in this dissertation. In the second study , a 
very simplified one-dimensional analysis through the thickness of the bolt hole 
was conducted to predict only the interlaminar normal stress due to combined 
thermal and bearing loading, neglecting all other stress components. The final 
study 5 3 was a two-dimensional study that showed in- plane stress concentrations 
at the bolt hole due to thermal loading. Hence, the study presented here on the 
thermal residual stresses in the vicinity of a bolt hole in a composite laminate is 
believed to be of much greater resolution then has previously been attained. In 
the present study, valid interlaminar stress concentrations at a bolt hole due to 
thermal loading and the manner in which the residual stresses change with a 
change in temperature, are investigated for the first time. 

5.2.1 Residual Stress Analysis Model 

The analysis of residual stresses was initially performed on only the 
section on the computational model given in Figure 5.2(b). The analyses were 
performed subsequently on the entire finite element model given in Figure 5.2(a), 
where no variations on residual stress components in the partial section were 
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observed between analysis models. This result shows that the residual stress 
concentrations for the current problem are confined to the region modeled in 
section b), and this section is sufficient for predicting residual stresses in the 
vicinity of the bolt hole for laminates of that size and larger. 

Boundary Conditions. For the residual stress analysis of the complete model, the 
boundary condition constraints included: 1) symmetry about z=0, where 
displacements in the z-direction were constrained, 2) symmetry about y=0, where 
displacements in the y-direction were constrained, and 3) displacements in the x- 
direction along the x = -1. inch boundary were constrained. 

Thermal Loading. The three operating temperature cases of 350°F, 70°F, and 
-200°F were applied as thermal loads to the analysis model with corresponding 
material properties, all having a reference stress free temperature of T S ft = 
460°F. A linear static analysis was performed for each thermal load case. 

5.2.2 Analysis Results 

Linear static thermal-structural analyses are performed to predict the 
thermal residual stresses in the vicinity of the bolt hole due to curing. At the bolt 
hole boundary, the in-plane residual stresses at T=70°F, predicted around the 
circumference of the bolt hole, are presented in Figure 5.3. For each stress 
component, the distribution presented is at the z-coordinate location through the 
thickness of the laminate where the maximum of that stress component occurred. 
Note that for o y , only the material matrix direction stresses were considered for 
the determination of the maximum in this study. This resulted in the maximum a y 
occurring in a resin layer. For the in-plane stresses, the most significant stress 
concentration due to the bolt hole occurred with a x , where a maximum of -32.ksi 
was observed at the bolt hole boundary. The cr x stress at the bolt hole boundary 
was twice the c x due to curing in the laminate without a bolt hole, where a value 
of -15. ksi was determined and presented in Table 3.2. 

The interlaminar stress distributions due to curing, determined at T=70°F, 
are presented in Figure 5.4 around the bolt hole boundary. Once again, for each 
stress component, the distribution is presented at the z-coordinate location 
through the thickness of the laminate where the maximum of that stress 
component occurred. As can be observed, very large interlaminar stresses occur 
at the bolt hole boundary, where each component has its peak at a different 
circumferential location. Table 5.1 presents the peak interlaminar stresses at the 
bolt hole boundary along with the peak interlaminar stresses at the straight free 
edge that were predicted in Chapter 3. As can be observed from Table 5.1, 
residual stresses due to a bolt hole boundary as compared to a straight free edge 
are dramatically larger. The o z stress concentration is over twice as large at the 
bolt hole boundary compared to the straight free edge. Where Ozx is zero for a 
straight free edge, a o zx of 16. ksi is predicted at the bolt hole boundary, and 
where a maximum positive o yz of 11.5 ksi is predicted at a straight free edge, a 
larger negative Cyz of -17.9 ksi is predicted at the bolt hole boundary. At the 
other operating temperatures studied, the circumferential stress distributions 
were similar with maximum values occurring at the same location. The only 
difference observed was in the magnitude of the stress components. 
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All peak stress components predicted in the laminate containing a bolt 
hole, at each operating temperature, are presented in Table 5.2. The variation of 
the maximum in-plane stress, o x , in the radial direction at 0 = 90 , from the bolt 
hole boundary, is presented in Figure 5.5, as predicted at the three operating 
temperatures being studied. Away from the bolt hole boundary, the stress 
concentration quickly reduces to the value predicted in a laminate without a bolt 
hole, where the values without a bolt hole are given in Table 3.2. Similarly, the 
variation of the maximum interlaminar stresses in the radial direction are 
presented in Figure 5.6, 5.7 and 5.8 for the o z , Oxz, and 0yz components, 
respectively. All interlaminar stresses quickly reduce to zero within a small 
distance from the bolt hole boundary. The behavior for the cure residual stresses 
at a bolt hole boundary is similar to the behavior of cure stresses in a laminate at 
a straight edge, where the interlaminar stresses also approach zero within one 
laminate thickness away from the free edge. 

5.3 Fastener Fit 

When the pin bearing configuration, as illustrated in Figure 1 .4, is exposed 
to a change in operating temperature, both the bolt hole diameter, depth and 
shape and the pin diameter may change due to either thermal expansion or 
contraction. The aforementioned pin-bearing configuration changes affecting the 
fastener fit have the potential to change the actual pin bearing surface area of the 
composite laminate. A change in the pin-bearing surface area can alternatively 
be caused by a change in the manufacturing clearance between the fastener and 
pin. The effect of clearance has been previously studied in the literature, and 
changing the pin-bearing surface contact area has been shown to affect the pin- 
bearing load capacity. 54 Additionally, pin bearing strength has been shown to be 
dependent on the bolt hole diameter size. 5 5 Hence, the study presented here is 
focused on determining the change in both the bolt hole diameter size, AD, and 
the bolt hole clearance, X, due to a change in temperature. 

The change in the diameter of the isotropic pin can easily be determined 
as a plane strain problem of a cylinder subject to a uniform change in 
temperature. The unrestrained change in pin diameter due to a change in 
temperature is computed as 

Ad = a p ATd (5-1) 

where exp is the coefficient of thermal expansion of the pin, AT is the change in 
temperature, and d is the original diameter of the pin. The pin properties used in 
the current study are presented in Chapter 4. 

The change in the geometry of the bolt hole of the composite laminate due 
to a change in temperature was predicted in the finite element analyses 
conducted for determining residual stresses in section 5.2. Analysis results 
revealed that for the [90,0] 8s IM7/PETI5 lay-up and the current operating 
temperatures considered, the changes in bolt hole geometry were extremely 
small. For the cure cycle from the cure temperature of T S ft = 460°F down to a 
room temperature of T = 70°F, a maximum change in the bolt hole diameter of 
only -2.82 x 10 -4 inch was predicted. The variation is only a 0.1% change in 
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diameter and is much smaller than the bolt hole tolerance of 0.001 inch used for 
machining the pin-bearing specimens. Additionally, a maximum change in 
laminate thickness of only - 8.19 x 10' 4 inch was predicted, which is only a 0.9% 
change and is much less then the variation in manufacturing thickness of the 
laminates, which were as large as 0.02 inch. The clearance, X, resulting from a 
change in temperature from a room temperature of 70°F to another operating 
temperature, is defined as 

X = AD - Ad (5.2) 

where AD is the maximum bolt hole diameter change computed from the 
differences between the finite element analysis results, and Ad is pin diameter 
change computed from Equation 5.1. 

The pin-bearing configuration changes affecting fastener fit, determined 
for the elevated and reduced operating temperatures being studied, are given in 
Table 5.3. Presented in Table 5.3 are the temperature change, AT, pin diameter 
change, Ad, laminate diameter change, AD, resulting clearance, X, and laminate 
thickness change, At. All diameter changes are extremely small relative to the 
dimensions of the pin-bearing problem. The net resulting clearances are also 
small, where for the elevated temperature of 350°F, a small decrease in 
clearance of -1.84 x 10 -4 is determined, and at the -200°F a small increase in 
clearance of 4.06 x 10 -4 is predicted. In relation to the effect of change in 
clearance on changes in pin bearing load capacity, a previous study 54 used a 
clearance of 0.01 and determined a decrease in load capacity of 12%. 
Therefore, the much smaller clearances predicted here due to temperature 
changes are assumed to have an insignificant effect on pin-bearing strength. 
Overall, changes in fastener fit due to the temperature changes of the current 
problem being studied are determined to be negligible and consequently will not 
be considered in the forthcoming pin bearing analysis models. 

5.4 Pin-Bearing Analysis 

A pin-bearing analysis model is developed where much of the detail is 
based on the findings of the previous studies presented in this dissertation. The 
effect of temperature changes on the state of stress at the bolt hole will be 
investigated for the pin-bearing load condition. The significant factors influencing 
the change in bearing strength with a change in temperature will be investigated. 

5.4.1 Load and Boundary Conditions 

Several methods of applying the pin-bearing load to the finite element 
model were considered to determine the optimum loading condition that could be 
modeled using NASTRAN. Many methods used in the literature were 
considered. The conditions to be discussed include: applied cosine load at the 
bolt hole boundary, fixed displacements at the bolt hole boundary, plate elements 
in the bolt hole, radial constraint at the bolt hole boundary, circumferential contact 
area, and nonlinear contact conditions including friction. 

The method of applying a cosine distributed load on the loaded half of the 
bolt hole boundary has been used in the literature on two-dimensional models of 
composite laminate. 5657 However, the method proved to be inaccurate for 
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composite problems with highly orthotropic properties. 5 8 The method was tested 
on the current problem and resulted in unrealistic stress concentrations at the 
interfaces between plies. Some authors have completely fixed the 
displacements on the loaded half of the bolt hole while applying a uniform force 
or displacement to the far end of the analysis model. 5 7,5 9 The method simulates 
a rigid pin with infinite friction, as in the case of a bonded pin. One study showed 
the method to result in unrealistic displacements for the bolt hole. The method 
was also tested on the current problem. At 9 = 90°, the intersection between the 
fixed displacement constraint and free bolt hole boundary, unrealistically high 
tensile stresses were observed. To account for pin elasticity, other authors have 
modeled the pin with bars and plate elements. 5 19 ' 511 However, for problems with 
zero clearance and high stiffness pins, the pin elasticity effect was negligible as 
predicted from both two- and three-dimensional nonlinear contact analyses. 

Since a high stiffness pin and zero clearance are the conditions of the current 
problem, the effect of pin elasticity should not be a significant factor and will 
consequently not be considered in the current pin-bearing analyses. 

The effect of friction between the pin and composite laminate has been 
studied using both two- and three-dimensional nonlinear contact analyses. 

As observed in the studies using a coefficient of friction of 0.2, the effect of 
friction has been to decrease the normal contact traction at 0 = 0 by about 25% 
while inducing a tangential contact traction. As 0 increases from 0° to 40° around 
the circumference of the bolt hole, the stress results approach the same stresses 
as were predicted in the case of no friction. Three studies 5 12,5 were found 
in the literature that performed three-dimensional nonlinear contact analyses of 
the pin bearing problem including friction. All the analyses used very crude 
meshes with a maximum of four ply layers modeled, and they all did not show 
any verification of the accuracy of the solution in the vicinity of the bolt hole 
boundary. The three-dimensional contact problem is highly nonlinear, requires 
significant computational resources, and may easily become unstable and 
produce erroneous results. 515 Since the current model contains an extremely 
large number of elements, modeling nonlinear three-dimensional contact is not 
practical. Also, since temperature effects on the pin-bearing problem are being 
studied for the first time in a three-dimensional high resolution model, a more 
reliable and efficient boundary condition was desired. 

The boundary condition of constraining only the radial displacements on 
the loaded half of the bolt hole was used to simulate a rigid frictionless pin. 
The load was then applied to the straight edge of the composite as a uniformly 
distributed load. This procedure was adopted for the current analysis as the 
most appropriate method of modeling the pin-bearing load. An initial analysis 
was performed where the radial constraint was applied from 0=0° to 0=90° on the 
bolt hole boundary. This condition resulted in erroneous tensile radial stresses at 
0=90°, which would not be the case of a laminate free to separate from the pin at 
this location. It was desired for the solution to only contain negative radial 
stresses in contact with the pin, which would be a realistic scenario of pin 
contact. The radial constraint was subsequently reduced by applying it to only 
the bolt hole boundary between 0=0° and 0=80°, simulating the case for a pin 
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having a slightly smaller diameter than the bolt hole, i.e., a pin with a clearance 
fit. This resulted in compressive radial stresses for pin contact between 0=0° and 
0=80°. 


Since the effect of friction was observed to be most significant only 
between 0=0° and 0=40° in a previous nonlinear contact analyses, the modeling 
of infinite friction was investigated here by constraining all displacements at the 
bolt hole boundary for the case of infinite friction between, 0 = 0° to 45° and in 
another case between 0 = 0° to 10°, while continuing to constrain only the radial 
displacement to 0 = 80°. Both cases resulted in unrealistic tensile stresses at the 
intersection between the fixed constraint and radial only constraint, which was 0 
= 45° in the first case and 0 = 10° in the second case. Hence, only radial 
constraints between 0=0° and 80° produced acceptable boundary stresses and 
these boundary conditions are used. 

In summary, the pin-bearing boundary and loading conditions used in the 
current analysis model are illustrated in Figure 5.9 and include 

1) radial displacement constraints applied to the bolt hole boundary 
between 0 = 0° and 0 = 80° 

2) y-displacement constraints on y=0 plane 

3) z-displacement constraints on z=0 plane 

4) uniform surface traction S g at x = -1. inch. 

The first boundary condition simulates the conditions of a rigid pin by 
constraining radial displacement at the bolt hole boundary and assumes a 
frictionless interface between the pin and composite laminate by leaving 
tangential and z-displacements unconstrained. The y- and z- model symmetry 
are prescribed with the second and third boundary conditions, and finally the load 
is applied in the negative x-direction by Sg. 

For the pin-bearing analyses, a pressure load equivalent to the pin- 
bearing yield stress, F bry , is applied for the evaluation of the state of stress in the 
linear range. Therefore, S g is determined as 


•Sg " 


FjryV 

w 


(5.3) 


where D = 0.25 inch is the bolt hole diameter, and W = 1 . inch is the finite 
element model laminate width. For each operating temperature analyzed, the 
corresponding average F bry as determined from experiment (Table 4.2) was used 

for determining the S g used in the analysis as presented in Table 5.4. 


5.4.2 Pin-Bearing Analysis Results 

Pin-bearing analyses were conducted to determine the state of stress at 
the three operating temperatures of 350°F, 70°F, and -200°F due to the pin- 
bearing load condition. For the room temperature operating condition, T= 70°F, 
in-plane stresses around the circumference of the bolt hole are displayed in 
Figure 5.10. The in-plane stress component plotted is at the z-coordinate 
location were the maximum occurred for that stress component. The interlaminar 
stress components are displayed in Figure 5.11, once again at the z-coordinate 
location were the maximum occurred for that stress component. 
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The distributions shown in Figures 5.10 and Figure 5.11 are very similar 
in shape to the distributions plotted in Ref. 5.16, where a three-dimensional pin- 
bearing analysis of a [90 4 ,0 4 ]s graphite/epoxy laminate was conducted using 
similar boundary conditions. The boundary condition differed from the current 
analysis in that the radial constraint was applied at the bolt hole between 0 = 0° 
and 0 = 90°. Both analyses had similar a x distributions with a compressive 
maximum at 0 = 0° and a tensile maximum at 0 = 90°. Where Ref. 5.16 shows a 
tensile o z peak at approximately 0=90°, a peak positive at 0=45°, and a peak 
negative o yz at approximately 0 = 22°, the current analysis model produced a 
peak o z at approximated 0 = 85° a peak axz at 50°, and, a peak Oyz at 27.5 . 
Overall, the similarity in appearance of the distributions gives validity to the 
current analysis model. 

The peak values for all stress components predicted from the analysis 
models at all three operating temperatures are given in Table 5.5 and the 
predicted bearing yield strain is presented in Table 5.6. Also presented in Table 
5.6 are the experimentally determined bearing yield strain, which will be 
discussed later in section 5.5. For all three operating temperatures, analysis 
results are plotted together for o x , a z , and Oyz in Figures 5.12-5.14, respectively. 
For each stress component, the distribution is presented at the z-coordinate 
location where the maximum value of that stress component occurred. The 
stress components plotted were chosen to show the similarity in the distributions 
between the three operating temperatures analyzed. The distributions were also 
similar for o y , and o xz , between the three operating temperatures, which were 
not plotted. It is important to note that the cause for variations in the state of 
stress are the variations in the lamina properties used for the analysis models at 
the different temperatures, in addition to the change in the magnitude of the load, 
S g , used in the analyses. Most importantly, variations in lamina properties with a 
change in temperature affect the state of stress at the bolt hole boundary and 
consequently should affect changes in the pin-bearing strength with a change in 
temperature. 

Additionally, to check the effect of modeling the full size specimen, an 
analysis was performed on an enlarged finite element model as can be seen in 
Figure 5.15. The full size model contained 80,310 elements, where the 
discretization in the 0.25 in by 0.25 in region adjacent to the bolt hole was not 
altered from the previous model given in Figure 5.2. A pin-bearing analysis was 
performed at the 70°F operating temperature, and the applied load was 
accordingly adjusted to S g = 6000. ksi. Except for the in-plane shear stress and 
interlaminar normal stress, the results showed less than 1% change in maximum 
stress magnitudes and no change in location. However, a 9% change in cr xy was 
observed were the maximum stress changed from -29. ksi to -32. ksi at 0 = 
22.5° and z=0.0818 in. Also, the maximum o z magnitude changed by 12% from 
18. ksi to 16. ksi, but the location or distribution was not altered. Although the 
changes in stress are significant in magnitude and indicate that the larger model 
results would be more appropriate for application of failure theories, the findings 
of the dissertation are not altered, and so the smaller model being used is 
concluded to be adequate for the current study. 
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5.4.3 Combined Analysis Results 

The state of stress due to the combined effect of pin-bearing load 
conditions and thermal residual cure stresses is obtained by superposition of 
solutions. The initial state of stress obtained from the analysis due to thermal 
cure loads at the operating temperature of interest is combined with the stress 
results obtained from the pin-bearing load conditions applied to the finite element 
model. 

At the room temperature operating condition of 70°F, the initial cure 
stresses are plotted along with the combined stress results around the 
circumference of the bolt hole for each stress component in Figures 5.16-5.21. 
Each stress component has been plotted at the z-coordinate location where the 
maximum of that stress component occurred. Unlike the in-plane stresses, all 
interlaminar stress components had the maximum occur at the same location as 
in the case of pin-bearing load only. The presence of the residual cure stresses 
only affected the magnitudes of the interlaminar stress component, and not the 
location of the maximum. Where the presence of cure stresses resulted in a 
lower maximum for the c x stress component as can be observed in Figure 5.16, 
this was not the case for the other stress components. For all other stress 
components, the magnitude of the peak stresses were intensified with the 
presence of the thermal residual cure stresses. For all interlaminar stresses 
plotted in Figures 5.19-5.21, the 0 location where the peak thermal cure stress 
occurred was near the same location where the peak pin-bearing stress 
occurred, resulting in intensified interlaminar stresses due to residual cure 
stresses. This effect is contrary to what was observed in Chapter 3 for the 
laminate, where the presence of residual stresses reduced the interlaminar 
stress magnitudes that were predicted under compressive loading. 

For the other operating temperatures of 350°F and -200°F, analysis 
results were similar, where changing the operating temperature did not affect 
either the location of the maximum or the appearance of the distribution. 
Therefore, their distributions were not plotted. Only the magnitude of the 
combined stresses varied with varying temperature. The maximum stresses 
predicted for each stress component and for each operating temperature are 
presented in Table 5.7. 

5.5 Comparisons with Experimental Results 

The computational results obtained in this chapter are evaluated through 
comparisons with experimental results. 

The magnitudes of the peak residual stress components are presented in 
Table 5.2. Through comparison with lamina strengths given in Table 2.1, their 
values reveal significant interlaminar stresses at all temperatures due to curing. 
The maximum o z exceeds the F^ strength predicted at all temperatures while o yz 
and Oxz exceed the Fi z strength at the operating temperatures of 70°F and 
-200°F. However, the strengths given in Table 2.1 only represent the average 
strength of the cross-section used in the property test and neglect stress 
concentrations in the test specimen. Note that the predicted maximum o z is 
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tensile, and hence a tensile modulus, E 33 T, would have been more appropriate in 
the elements where o z were positive. Since E 22 T is lower then the E 22 C modulus 
used in the analysis, and E 33 =E 22 is assumed in the model, lower stresses would 
be expected if a tensile modulus were utilized. However, a qualitative 
perspective can be gained from the analysis results presented. 

The magnitudes of the maximum stress components are presented in 
Table 5.5. Through comparison with lamina strengths given in Table 2.1, their 
values reveal significant stresses for all stress-components due to the pin- 
bearing load condition. For pin-bearing loading only, the in-plane shear stress 
and all interlaminar stresses exceed the strength values. Consequently, the case 
of pin-bearing loading only over-predicts stresses at the bolt hole in the 
composite laminate. 

The bearing yield strain predicted from the pin-bearing analysis is 
compared to the measured values from the experiment in Table 5.6. Large 
discrepancies are observed. Unfortunately, the previous studies that used the 
same method for measuring bearing strain did not attempt to verify their results 
with an analysis as was attempted here . 5 17,518 There is a possibility here that 
the measured values are erroneous and that the extensometer used did not 
measure the bolt hole elongation accurately. Although the measured values are 
questionable, the bearing strain measurement is still useful and necessary for 
indication of the yield and ultimate strength values obtained from the experiment. 
Also, an observation can be made that both the measured and predicted bearing 
strains were lower at room temperature than at the elevated or reduced operating 
temperatures. 

The maximum stress failure criterion is used here to gain perspective on 
the severity of the stresses predicted from the analysis. The maximum stress 
criterion states that each stress component must be less than the respective 
strength, otherwise, failure is said to have occurred . 519 A failure ratio is 
computed as 

F.R. = — (5.4) 

F, 

where ay is the maximum stress for the ij component and Fy is the strength of the 
corresponding stress component. Then, a failure ratio of F .R. > 1 indicates 
failure. The failure ratios computed from the maximum stresses are given in 
Table 5.8. All in-plane shear stresses and interlaminar stresses are well above 
the failure ratio, indicating failure had occurred before the applied loads were 
reached. Consequently, either the strengths have been under-predicted from the 
experiments or the model has over-predicted stresses. 

With respect to the strength data, the interlaminar strengths given in Table 
2.1 do not account for stress concentrations during the experiment and may 
actually under-predict the strength in the absence of stress concentrations. Also, 
since resin layers are present even in unidirectional laminates, residual stresses 
can also be present and affect the lamina strength data, where a correction could 
be made to the lamina strength data to account for residual stresses. 

With respect to the stresses predicted by the model, several reasons can 
be offered for over-prediction. First, as mentioned earlier in evaluation of 
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residual stresses, the tensile modulus E 22T is lower than the compressive 
modulus E 2 2 c at all operating temperatures (see Table 2.1). Consequently, 
appropriate use of the tensile E 22T , in elements where a tensile stress was 
predicted should lead to lower tensile stresses. The high failure ratio for a xy at the 
350°F operating temperature could be lessened by appropriately modeling the 
nonlinear material behavior, as obviously would be necessary in this case. 
Finally, the effect of friction at the bolt hole boundary, as was neglected here 
could result in lower boundary stresses. 

In conceding to over-prediction, however, additional observations can be 
made on the predicted state of stress. With the failure ratios being near the 
same value if not greater with increasing operating temperature, and with the 
lower average bearing yield strengths values being used as loading with the 
increasing operating temperatures, the analyses reveal that increasing 
temperature reduces the pin-bearing yield strength. This observation 
corresponds to the experimental results of decreasing strength with increasing 
temperature. 

The distribution of the maximum interlaminar stresses through the 
thickness of the laminate, at the bolt hole boundary and 0 location where the 
maximum occurred, are presented in Figures 5.22-5.24. Plotted are only the 
distributions for the case of T=70°F, since similar distributions were observed at 
the other operating temperatures. As seen in Figure 5.22, the maximum 
interlaminar normal stress occurs in the lower 90° ply at the interface with the 
lower resin layer. The peak layer stress is then observed to diminish as the outer 
plies are approached. In contrast, the interlaminar shear stress maximum, Oxz, 
occurs in the outermost resin layer interface with the 0° ply, and the maximum Oyz 
occurs in the second resin layer from the top interface with the outermost 0° ply, 
where both interlaminar stress layer peaks lessen as the center of the laminate is 
approached. In comparison with the failed specimen micrograph (see Figure 
4.8), the maximum interlaminar shear stress locations could be used to explain 
the severe damage observed in the specimens in the outermost two resin layers. 
Consequently, failure of the composite may initiate from excessive interlaminar 
shear stresses. 

To aid in understanding the contribution of thermal residual stresses on 
the maximum stress result, the magnitude of the cure stresses at the location 
where the maximum combined stress occurred are presented in Table 5.9. As 
can be observed from the table, significant changes in residual stress occur 
between T=350°F and T=70°F with a AT=280°F. The changes in residual stress 
are significant in comparison to the change in bearing strength, where a A F bry of 
only -6ksi was determined from experiment. Also, since the residual stresses are 
lower at T=350°F than at T=70°F, the effect of thermal residual stress is 
concluded to lessen the degradation in pin-bearing strength than would be 
expected in the absence of cure stresses. With the changes in residual stresses 
being more negligible between 70°F and -200°F, the change in pin-bearing 
strength to the reduced temperature appears to be affected by only the change in 
the lamina material properties. Here, a greater change in F bry of +15 ksi was 
measured for AT= -270°F. Perhaps the change in strength is greater since there 
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are no changes in thermal residual stresses lessening the change in strength as 
there was between 70°F and 350°F. 

The presence of residual stresses has been shown previously to intensify 
the magnitude of the maximum in-plane shear and interlaminar stresses obtained 
under pin-bearing loads for all temperatures. Consequently, residual cure 
stresses appear to contribute to the reason why the bearing strengths in general 
are lower than the compressive strength of the composite laminate. This is true 
for all lay-ups tested, where all laminate strengths presented in Table 2.4 are 
greater than the corresponding average ultimate bearing strengths presented in 
Tables 4.2, 4.4, 4.5, and 4.6. 

Most importantly, from the results of this analysis, in addition to a change 
in lamina material properties, changes in the residual cure stresses with a 
change in temperature have the potential to affect the change in pin-bearing 
strength with temperature change. Hence, both AM(A1) and AS(AT) are factors 
that cannot be ignored in an investigation on the change in bearing strength with 
a change in temperature, AF br (A7). 
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Table 5.1 Peak Interlaminar Residual Stresses Predicl 

ted at T=70°F. 

Stress Component 



o z (ksi) 

7. 

16.7 

Oxz (ksi) 

0 . 

16.0 

o yz (ksi) 

11.5 

-17.9 


Table 5.2 Peak Residual Stresses in Laminate Containing the Bolt Hole- 



Temperature 

Stress Component (ksi) 

350°F 

70°F 

-200°F 

Ox 

-6.9 

31.9 

31.7 

5* 

5.7 

21.9 

23.7 

tfxy 

-0.86 

-4.5 

-6.2 

O z 

3.5 

16.7 

16.1 

Oxz 

3.2 

16. 

16.7 

Oyz 

-3.6 

-17.9 

-18.7 


Table 5.3 Changes in Fastener Fit Due to Temperature Changes. 


Temperature 

AT 


AD 

X 

mu 

(°F) 



(inches) 

(inches) 


350 

280 

3.92 x 1 0' 4 

1.94 x 10 4 

-1.98 x 1 0' 4 

5.8 x 10' 4 

-200 

-270 

-3.78 x IQ- 4 

0.14 x IQ- 4 

3.92 x 1 0‘ 4 

0 . 


Ta ble 5.4 The Pin-bearing Load, Sq, for each Operating Tempera 




■u@m 

mMSENKm 



-200 

63. 

-15.75 

70 

48. 

-12.00 

350 

42. 

-10.50 


ure 


Table 5.5 Maximum Stresses due to Pin-Bearing Loads. 



Temperature 

Stress (ksi) 

350°F 

70°F 

-200°F 

o x 

214. 

206. 

265. 

Oy 

-30. 

-31. 

-21. 

&xy 

-26. 

-29. 

-41. 

o z 

10. 

18. 

14. 

Oxz 

24. 

23. 

31. 

GyZ 

35. 

35. 

44. 
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Table 5.6 Bolt Hole Bearing Yield Strain. 


Temperature 

- ■ '- r '-t 

Experiment 

Analysis 

(°F) 

Ebrv (j-iin/in) 

Ebrv (H'n/in) 

350 

30300. 

11400. ' 

^ 70 

21700. 

10920. 

-200 

23500. 

13680. 


Table 5.7 Maximum Stresses from Combined Solution - 
Pin-Bearing Loads and Thermal Residual Cure Stresses. 



Temperature 

Stress (ksi) 

350°F 

70°F 

-200 F 


211. 

189. 

246. 

Oy 

22. 

45. 

44. 

Oxy 

-29. 

-35. 

-49. 

o z 

12. 

27. 

24. 

&XZ 

26. 

34. 

41. 

Gyz 

37. 

47. 

55. 



Temperature 

F.R. 

350°F 

70°F 

-200 F 

o x 

.69 

.62 

.73 

Oy 

.82 

1.22 

1.07 

Oxy 

4.39 

2.40 

2.75 

o z 

3.70 

6.00 

4.62 

Oxz 

2.17 

1.97 

1.92 

Gyz 

3.08 

2.72 

2.58 
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Table 5.9 Cure Stresses at Combined Solution Maxima. 



Temperature 

Stress (ksi) 

350°F 

70°F 

-200°F 

Ox 

-1. 

-17. 

-18. 

Oy 

3. 

15. 

15. 

Oxy 

-1. 

-5. 

-5. 

O z 

2. 

9. 

9. 

Oxz 

2. 

10. 

11. 

Oyz 

2. 

12. 

11. 



Figure 5.1. Pin-bearing composite laminate computational model geometry. 
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0.25 in 


a) complete finite element model shown in x-y plane. 



0.125 in 


b) Section of finite element model showing discretization 
at z=0 in x-y plane. 

Figure 5.2. Finite element model discretization in x-y plane. 
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0, degrees 


Figure 5.3. Variation of in-plane stresses around circumference of bolt hole 
due to curing, T = 70°F. 
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Figure 5.4. V 
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Figure 5.6. Radial distribution of o z at 6 = 0 degrees due to curing, at 
T = -200°F, 70°F, 350°F. 
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Figure 5.7. Radial distribution of Gxz at 9 = 52.5 degrees due to curing, at 
T = -200°F, 70°F, 350°F. 
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Figure 5.8. Radial distribution of Oyz at 0 = 37.5 degrees due to curing, at 
T = -200°F, 70°F, 350°F. 
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Figure 5.10. Variation of in-plane stresses around circumference of bolt hole 
due to pin-bearing load condition at T=70°F. 
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Figure 5.1 1 . Variation of interlaminar stresses around circumference of bolt hole 
due to pin-bearing load condition at T=70°F. 
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Figure 5.15. Finite element model for model size evaluation. 
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Figure 5.16. o x around circumference of bolt hole at z - 0.0821 inch. 



Figure 5.17. o y around circumference of bolt hole at z - 0.0083 inch. 



Figure 5.18. o xy around circumference of bolt hole at z - 0.0802 inch. 
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Figure 5.22. Combined solution o z distribution at bolt hole boundary and 
9 = 85° for T=70°F. 
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Figure 5.23. Combined solution Oxz distribution at bolt hole boundary and 
0 = 52.5° for T=70°F. 
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Figure 5.24. Combined solution Oyz distribution at bolt hole boundary and 
9 = 27.5° for T=70°F. 




Chapter 6 
Conclusions 
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An overview of the dissertation research accomplishments and findings is 
presented within this chapter. First, a summary reviewing the major results is 
presented. Then, specific conclusions from the research are presented. Finally, 
recommendations for future work are given. 

6.1 Summary 

Advanced technology vehicles require numerous composite structural panels 
to be mechanically joined and to operate at temperature extremes. Experimental 
studies have shown that elevated temperatures reduce the strength of polymer 
joints to an extent not explained by conventional understanding of joint behavior. 
A review of previous research did not identify a comprehensive study showing 
the factors contributing to the strength change. Consequently, the goal of this 
dissertation is to develop a fundamental understanding on the mechanisms that 
cause strength changes with a change in temperature. A fundamental 
understanding is necessary to design an efficient composite joint able to 
withstand a change in operating temperature without the risk of premature failure. 
A hypothesis was developed relating the change in pin-bearing strength to three 
basic changes in a specific joint subject to temperature change. The proposed 
factors contributing to a strength change with a change in temperature include: a 
change in the state of residual thermal stress, a change in the material properties 
of the constituent materials, and a change in the fastener fit. The specific 
objectives are then given with the goal of the research being to investigate the 
significant factors contributing to a strength change with a change in 
temperature. The objectives include: characterizing the material system to model 
structural response accurately, conducting pin-bearing experiments under 
thermo-mechanical loads, and developing computational models to study the 
effect of temperature changes on the state of stress. 

The material system characterization to model the IM7/PETI5 material 
system utilized in this study is presented. The temperature-dependent lamina 
engineering constants were determined between temperatures of 495°F and 
-200°F. The glass transition temperature of the IM7/PETI5 was determined to 
be 460°F, and this temperature was assumed to be the stress free temperature 
below which curing stresses develop in the laminate. Lamina strength properties 
were determined at operating temperatures of -200°F, 70°F, and 350°F. 
Nonlinear stress-strain behavior was presented for the in-plane shear modulus 
and the matrix-direction compressive modulus. Laminate strength properties for 
several lay-ups of interest were also determined for comparisons with 
computational results and pin-bearing test data. 

As a preliminary step in developing a computational model of the pin- 
bearing problem, analysis results for the composite laminate were studied. 
Classical lamination theory, as pertained to the current problem of interest, was 
introduced, and the derivation of the classical lamination theory expression for 
computing in-plane stresses due to both mechanical in-plane loads and a uniform 
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temperature change was presented. Key background information on interlaminar 
stresses and residual cure stresses revealed the need for a three-dimensional 
analysis of the composite laminate. Accordingly, a three-dimensional finite 
element model of a composite laminate was developed using NASTRAN and 
validated with an analysis found in the literature. The model was used to 
investigate laminate attributes that significantly affect the state of stress in the 
laminate and hence need to be incorporated in the pin-bearing analysis. The 
need for modeling the resin layers between lamina and residual stresses in the 
laminate was established. For the operating temperatures under investigation, 
using constant material properties determined at the operating temperature of 
interest in a linear analysis provided an accurate stress solution, as opposed to 
modeling temperature-dependent properties and conducting nonlinear analyses. 
Also, for the [90,0] 8s lay-up, analysis results revealed the in-plane shear stresses 
and interlaminar shear stresses, xi 3l to be virtually zero, eliminating the need for 
nonlinear stress-strain characterization of the laminate model. 

A combined solution was produced where the residual stresses due to 
curing were combined with the solution of the composite laminate subject to a 
compressive load. The interlaminar stresses present in the laminate due to 
curing were found to be reduced with the application of the mechanical 
compressive load. Therefore, the cause of failure in a composite laminate 
subject to compressive load is most likely not due to excessive interlaminar 
stresses. Results of the combined analysis also revealed that a change in 
laminate strength with a change in temperature is governed by the change in the 
fiber-direction compressive lamina strength, Fn c , with a change in temperature. 
The possibility that the residual stresses present in the laminate may lessen the 
change in strength was explored. Moreover, a three-dimensional analysis was 
concluded to be necessary to determine the role of interlaminar stresses in the 
pin-bearing problem. 

Pin-bearing experiments were conducted to produce pin-bearing strength 
data for developing understanding on the pin bearing behavior of the IM7/PETI5 
material system. The tests were conducted at the three operating temperatures 
of interest, 350°F, 70°F, and -200°F. As expected based on previous studies, for 
all lay-ups that failed in the bearing mode, strength degradation was observed 
with increasing temperature. Only the [90]e s lay-up, which failed in the tension 
mode, did not show any significant strength change with changing temperature. 
In comparing laminate and pin-bearing strength changes with temperature 
change, the extent of the strength change differs between the laminate and pin- 
bearing strength and also varies with lay-up change. Also, most pin-bearing 
strengths were significantly lower than laminate strengths. Since the 
computational study was limited to the [90,0] 8s lay-up, this lay-up was scrutinized 
more thoroughly than the others. For all lay-ups, the ultimate bearing strengths 
and failure modes were presented. For the [90,0] 8s lay-up, the bearing yield 
strength, bearing yield strain, and bearing ultimate strain were also computed 
and presented for all specimens tested. Micrographs taken on failed specimens 
cut along the bearing plane revealed severe damage in the outermost two plies. 
Also emanating from that region were transverse shear cracks. This observation 
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raised the possibility that ultimate failure could be due to the interlaminar shear 
stresses in the region of severe damage. 

The computational study of the pin-bearing problem was finally presented. 
A three-dimensional finite element model of the composite laminate containing 
the bolt hole was developed. The model was highly discretized at the bolt hole 
boundary and contained a total of 63,498 elements. The residual stresses 
present in the laminate due to curing from the assumed stress free temperature 
of 460°F down to the operating temperature were studied. The free edge 
stresses at the bolt hole boundary were found to be dramatically different than in 
the case of a straight free edge. All interlaminar stress components intensified 
with extremely large stress concentrations at the hole boundary. As with the 
straight free edge, the interlaminar stresses quickly reduced to zero within one 
laminate thickness away from the free edge boundary. The residual stress 
distributions were similar for all three operating temperatures studied. Where the 
magnitudes of the maximum stresses of each component were lowest at 350°F, 
there was not much difference in residual stresses between the 70°F and -200°F 
operating temperatures. The effect of fastener fit changes with a change in 
temperature was also studied and found to be negligible for the current pin- 
bearing problem. 

The finite element model was analyzed subject to pin-bearing loads. 
Assuming a rigid frictionless pin, the pin-bearing load was modeled assuming 
radial displacement constraints at the bolt hole boundary between 0 = 0° and 80°, 
and a uniform negative pressure load to the straight end of the model at x = -1. 
inch. For each temperature analyzed, the load applied was equivalent to the 
average pin-bearing yield strength determined from the experiments of Chapter 
4. Analysis results revealed extremely high stresses for all stress components at 
different locations around the bolt hole boundary. Distributions were similar for 
all three operating temperatures analyzed, where variations in only the 
magnitude of the maximum stresses were observed. Variations in material 
properties were obviously concluded to affect the magnitude of stresses 
significantly at the bolt hole boundary. 

Finally, a combined solution was obtained that included the thermal 
residual cure stresses in the state of stress for the pin-bearing problem. The 
presence of cure stresses in the laminate intensified maximum pin-bearing 
stresses at the bolt hole boundary for all operating temperatures studied. 
Consequently, residual stresses could possibly contribute to the reason why pin- 
bearing strengths are lower than compressive strengths of the composite 
laminate. The observation was also made that the presence of residual stresses 
may lessen the extent of changes in pin-bearing strength with a change in 
temperature. Failure ratios much greater than one were predicted for many of 
the stress components, indicating the possibility that the finite element model did 
not include enough detail to accurately predict stresses. However, the 
distributions predicted peak interlaminar stresses in the outermost two resin 
layers of the composite laminate at the bolt hole. The micrograph of the failed 
specimens also revealed severe damage in the outmost two layers, indicating the 
possibility of failure initiating in that region. 
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6.2 Conclusions 

The following conclusions were reached based on the experimental and 

computational studies investigating the effect of temperature changes on 

changes in pin-bearing strength of an IM7/PETI5 composite laminate: 

1. Changes in both residual stresses and material properties with a change in 
temperature are factors influencing the changes in pin-bearing strength with 
temperature changes as proposed in the hypothesis. 

2. The effect of changes in fastener fit with a change in temperature are 
negligible on pin-bearing strength changes and can be neglected from the 
hypothesis for the pin-bearing problem studied in this dissertation. 

3. A three-dimensional analysis model including the resin layers between plies 
and thermal residual stresses due to curing is required to investigate 
accurately the state of stress in the pin-bearing problem. 

4. The presence of residual stresses due to cure have the potential to reduce 
the overall pin-bearing strength of a composite laminate. 

5. Residual stresses due to curing have the positive affect of lessening the 
extent to which bearing strength decreases with a temperature increase from 
70°F to 350°F for the pin-bearing problem studied. 

6. Residual stresses do not significantly change when temperature decreases 
from 70°F to -200°F and material property changes are the dominant factor 
influencing the increase in pin-bearing strength for the pin-bearing problem 
studied. 

7. Analysis results reveal a decrease in the pin-bearing strength with increasing 
temperature, which qualitatively corresponds to the experimental pin-bearing 
strength reductions with increasing temperature. 

8. The finite element model did not include sufficient detail to predict stresses 
accurately for the pin-bearing problem. 

9. The reduced strength and stiffness properties of the polymer matrix material 
with increasing temperature drives pin-bearing strength degradation with 
increasing temperature and the residual thermal stresses due to curing only 
influence the extent of the strength degradation. The extent of the strength 
reduction with increasing temperature can only be lessened by the presence 
of the thermal residual stresses, which decrease with increasing 
temperature. 

10. Until an accurate pin-bearing computational model is developed to predict 
strengths accurately, an IM7/PETI5 (or similar advanced material) joint 
design, with a defined lay-up, needs to be tested at the desired operating 
temperature to determine its pin-bearing strength at that temperature. The 
upper bound operating temperature will determine the joints minimal pin- 
bearing load capacity. 


6.3 Recommendation 

The following recommendations for future studies are based on the findings 
of this dissertation. 
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A more accurate finite element model is needed to predict stresses for the 
evaluation of pin-bearing behavior including thermal residual stresses. First, the 
effect of friction on the magnitude of the stresses at the bolt hole boundary needs 
to be investigated. A three-dimensional nonlinear contact analysis capability is 
recommended for modeling the effect of friction. Secondly, due to the large 
variation of key stresses between very high negative and positive values, both 
the tensile and compressive moduli should be incorporated in the analysis. A 
non-linear iterative scheme, where the determination on whether the stress was 
positive or negative for an element is made with increasing load, so the 
appropriate modulus can be determined and utilized for that element, is 
recommended. Finally, due to the large shear stresses predicted for the pin- 
bearing problem, modeling of the nonlinear stress-strain behavior is 
recommended. The non-linear moduli need to be modeled independently, where 
for the IM7/PETI5 material, the fiber direction properties are linear, and the matrix 
and shear properties have differing nonlinear profiles. 

Due to the shortcomings of NASTRAN to include any of the aforementioned 
detail, an alternative analysis tool with the ability to model three-dimensional 
contact, tension and compression moduli simultaneously, and nonlinear moduli 
independently, needs to be identified, or if necessary, developed. 

With large discrepancies observed between the measured and computed 
bolt hole elongation, evaluation and if necessary modification of the 
measurement technique used in this dissertation is recommended. 



Appendix 

Material Property Testing 
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The tests conducted to determine the moduli and strength properties are 
outlined in Table 1 . The table shows the type of test followed by the test method 
specimen size, the number of plies in the laminates tested and finally the 
material properties obtained from the test. The ASTM test methods are pu 
by the American Society of Testing and Materials (ASTM) in their Annual book of 
ASTM Standards. The lamina property testing was conducted under gran 
University of Wyoming (exluding Fn c ), while the laminate property testing and the 
F 11c property tests were conducted at NASA Langley Research Center (LaRC). 

The CLC method uses the Combined Load Compression(CLC) fixture 
developed at the University of Wyoming. 1 The CLC method was chosen Jl® , | ® 
over the ASTM standards for compression testing of composites, since samp 
[90,0 ] 8s specimens tested using the CLC method yielded hig er v i 
compressive strengths. In addition, there was increased simplicity of conducting 
the test using the CLC fixture. While using the CLC fixture in testing at LaRC for 
the File strength and the laminate compression strengths, it was discovered that 
additional clamp-up on the fixture was needed to prevent P re m 3 ture failure on 
the loaded ends of the specimen for the 350 F tests. An increase in tor 9 ue 
in-lbs (from the recommended 25. in-lbs given in the user instructions) was found 
to be sufficient for failure to occur in the gage section. It was also found that 
using a torque too high, i.e. greater then 25 m-lb for the -200 F tests, result 
premature failure occurring in the gage section adjacent to the fixture. 

The property testing was conducted at the three operating temperatures of 
— 200°F 70°F and 350°F and at an upper temperature of 495 F, which was 
preliminarily determined to be the glass transition temperature T g of the 
material. Later, the glass transition temperature, T g , of the n ^ rre !l t ' M7/ ^ E . t 
material was determined, by Wallace Vaughn , to be 460 F The test data 
presented in Chapter 2 are typically the average values of four replicate 
specimen When specimens were determined to have failed prematurely, 
data was not used in catenating the average value presented The coeffoents 
of thermal expansion presented in Table 2.1 were determined at NASA Langley 
Research Center by Wally Vaughn and Craig Ohlhorst. 
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1 Aerospace Technologist, Enviromental Interactions Branch, NASA Langley Research Center, Hampton, 
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Table 1 Material Property Test Information for Material Characterization 


Lamina Properties 


Type 

Method 

Specimen 
Size (inxin) 

# of plies 

Properties 

Axial 

Compression 

CLC** 

0.5 x 5.5 

20 

£nc, fiic 

Transverse 

Compression 

CLC** 

0.75 x 5.5 

32 

flic, f 11C.V23 

Axial 

Tension 

ASTM 

D3039* 

0.5x 11. 

8 

f22T, f22T,Vi2 

T ransverse 
Tension 

ASTM 

D3039* 

0.75x11. 

20 

f 1 1T, f 11T 

Shear 

ASTM 

D5379* 

0.5 x 3. 

20 

f 12, f 12 

ILSS + 

ASTM 

D2344* 

0.5 x 1. 

32 

fxz 

ILNS ++ 

i 

Modified 
ASTM C 297* 

1. x 1. 

32 

fzz 

Laminate Properties 

Compression 

CLC** 

0.5 x 5.5 

32 

£xxc, Fxxc 

Tension 

ASTM 

D3039* 

0.5x 11. 

32 

fxXT, fxXT 


ASTM D# — American Society for Testing and Materials Designation Number 
CLC - Combined Load Compression 
+ ILSS - Interlaminar Shear Strength 
++ ILNS - Interlaminar Normal Strength 




